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ABSTRACT 


Heat  transfer  measurements  were  made  for  filmwise  condensation 
of  R-113,  steam  and  ethylene  glycol  on  three  sets  of  finned  tubes 
which  differed  only  in  root  diameter.  The  fin  root  diameters  were 
12.7  mm.  19.05  mm  and  25.0  mm.  A  comparison  of  the  enhancement 
ratios  (based  on  constant  vapor- side  temperature  drop)  revealed  that, 
within  the  range  of  diameter  tested,  the  effect  of  root  diameter  was 
small.  Results  indicated  that  two  or  more  trends  may  exist  while 
increasing  root  diameter. 

With  the  exception  of  the  small-diameter  tubes  with  steam,  a 
change  in  root  diameter  did  not  effect  the  optimum  fin  spacing  for 
each  fluid  tested.  The  optimum  fin  spacing  for  the  small-diameter 
tubes  with  steam  was  approximately  2.0  mm,  while  the  medium-  and 
large- diameter  tubes  had  an  optimum  fin  spacing  of  1.5  mm.  The  opti¬ 
mum  fin  spacing  for  R-113  and  ethylene  glycol  was  found  to  be 
0.5  mm  and  1.0  mm,  respectively. 

A  comparison  between  the  outside  heat-transfer  coefficients  of 
the  medium-diameter  tubes  for  R-113  and  steam  with  past  NPS  inves¬ 
tigation  showed  an  agreement  within  ±3  percent  and  ±10  percent, 
respectively. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Heat  exchangers  play  an  intricate  role  in  the  design,  building  and 
mission  capabilities  of  today’s  naval  vessels.  Main  engine  and  auxiliary 
condensers  ensure  the  safe  operation  of  the  main  engines  and  auxiliary 
equipment  as  well  as  the  generation  of  electrical  power.  Air-condi¬ 
tioning  condensers  help  to  provide  the  climate  control  for  ship’s 
company  and  sensitive  equipment.  The  development  of  high-powered 
shipboard  sensors,  computers,  and  weaponry  has  stretched  present 
shipboard  air-conditioning  units  to  capacity.  As  a  result,  the  modern¬ 
ization  of  naval  vessels  has  required  several  tons  of  additional  cooling 
capacity.  Of  course,  this  addition  means  that  more  space  and  weight 
must  be  allocated  to  this  equipment.  The  Naval  Postgraduate  School 
(NPS),  in  association  with  David  Taylor  Research  Center  and  with  sup¬ 
port  from  the  National  Science  Foundation,  is  presently  involved  in 
research  aimed  at  reducing  the  size  and  weight  of  condensers.  These 
efforts  can  reduce  ship’s  tonnage,  save  fuel,  and  provide  space  for 
other  uses. 

B.  CONDENSATION 

There  are  two  modes  of  condensation:  dropwise  and  fllmwise.  The 
random  formation  and  departure  of  discrete  droplets  on  the  con¬ 
densing  surface  is  known  as  dropwise  condensation.  This  mode  is  very 
efficient  because  the  hot  vapor  is  able  to  come  in  closer  contact  with 
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the  cooler  condenser  surface  because  of  the  presence  of  very  small 
droplets.  Unfortunately,  sustained  dropwise  condensation  has  not  been 
practical  despite  decades  of  research.  Filmwise  condensation  is  a  pro¬ 
cess  in  which  the  entire  surface  is  covered  with  a  thin,  continuous 
film  of  condensate.  The  film  creates  a  larger  resistance  to  the  heat 
transfer  and  therefore  this  process  is  less  efficient.  However,  in 
designing  and  testing  a  condenser,  filmwise  condensation  is  desired 
as  this  is  the  most  likely  mode  to  occur.  The  high  efficiency  and  ran¬ 
domness  of  dropwise  condensation  are  not  predictable  and  could  lead 
to  over-speculation  of  the  condenser’s  capabilities.  Designing  for  film- 
wise  condensation  is  conservative  and  leaves  a  margin  for  error 
[Ref.  1J. 

C.  BASIC  HEAT-TRANSFER  EQUATION 

The  basic  heat-transfer  equation  used  to  describe  the  relationship 
between  temperature  difference  and  heat-transfer  rate  in  a  heat 
exchanger  is: 


Q  =  UA(LMTD) 


(1.1) 


where 

Q  =  heat- transfer  rate, 

U  =  overall  heat-transfer  coefficient, 

A  =  surface  area  for  heat  transfer  consistent  with  U,  and 
LMTD  =  log-mean-temperature  difference. 
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In  a  condenser,  the  value  for  the  log-mean-temperature  difference 
(LMTD)  is  dependent  on  the  characteristics  of  the  vapor/ condensate 
and  cooling  water  flow.  The  limits  of  this  parameter  are  set  primarily 
by  the  saturation  temperature  of  the  vapor  and  the  temperature  of  the 
cooling  water  sink. 


LMTD  = 


(T.-T,) 

In 

^sat  ”T, 

Lt^-tJ 

(1.2) 


where 

Ti  =  cooling  water  inlet  temperature. 

T2  =  cooling  water  outlet  temperature,  and 
Tsat  =  vapor  saturation  temperature. 

The  overall  heat-transfer  coefficient  (U)  is  dependent  on  cleanli¬ 
ness,  construction,  and  composition  of  the  condenser  tubes,  as  well  as 
the  fluid  properties  of  the  condensate  and  the  cooling  water.  It  is 
inversely  proportional  to  the  resistance  to  heat  transfer.  The  area  (A), 
consistent  with  the  coefficient,  is  usually  the  outside  area  of  the  tubes. 
This  area  can  be  increased  by  the  addition  of  fins.  However,  as  pointed 
out  by  Wanniarachchi,  et  al.  [Ref.  2]  and  others  [Ref.  3,4),  the  enhance¬ 
ment  to  heat  transfer  observed  after  adding  fins  is  caused  not  only  by 
the  increased  area  but  also  by  the  thinning  of  condensate  film  due  to 
surface  tension  forces.  The  change  in  film  thickness  affects  the  resis¬ 
tance  to  heat  transfer  and  therefore  affects  the  overall  heat-transfer 
coefficient.  On  the  other  hand,  surface  tension  forces  cause  a  thick 
condensate  film  to  be  present  between  the  fins  over  the  lower  portion 
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of  the  tube.  This  tends  to  create  a  poorer  performance  on  the  lower 
portion  of  the  tube.  Despite  this,  finned  tubes  have  been  demonstrated 
to  provide  considerable  enhancement  over  smooth  tubes  for  various 
fluids,  including  steam  [Refs.  2,  3.  5).  By  considering  the  LMTD  to  be 
almost  constant.  Equation  1 . 1  clearly  states  that  the  optimum  value  for 
heat  transfer  occurs  at  the  optimum  value  of  the  product  of  the  surface 
area  and  the  overall  heat-transfer  coefficient. 

D.  RESEARCH  AT  THE  NAVAL  POSTGRADUATE  SCHOOL 

A  systematic  study  to  determine  the  effects  of  fin  parameters, 
such  as  fin  spacing,  height,  and  shape,  tube  diameter  and  vapor  veloc¬ 
ity  on  condensation  heat-transfer  enhancement,  began  in  June  1982 
with  the  building  of  a  test  apparatus.  Krohn  [Ref.  6]  designed  and  built 
the  experimental  apparatus  to  study  enhanced  condensation  heat 
transfer  on  horizontal  tubes.  His  apparatus,  as  modified,  is  described 
in  Chapter  III  of  this  thesis.  Graber  [Ref.  7)  completed  the  instrumen¬ 
tation  and  calibration  for  the  apparatus  and  took  preliminary  data  on  a 
smooth  tube.  The  original  smooth  tube  had  thermocouples  embedded 
in  the  tube  wall  to  directly  measure  wall  temperature.  This  tube  was 
used  to  generate  a  correlation  for  the  inside  heat-transfer  coefficient 
which,  in  turn,  was  used  to  infer  the  outside  heat-transfer  coefficient 
by  subtracting  the  inside  and  wall  resistances  from  the  overall  resis¬ 
tance  (Equation  4.17).  In  December  1983,  Poole  [Ref.  8[  demonstrated 
the  use  of  the  modified  Wilson  plot  technique  that  compared  favorably 
with  the  direct  method  involving  wall  temperature  measurements. 
The  modified  Wilson  plot  technique  is  described  in  Chapter  IV  of  this 
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thesis.  Poole  experienced  difficulties  in  obtaining  complete  filmwise 
condensation  on  a  copper  tube  with  steam  as  a  working  fluid.  He  had 
partial  success  in  obtaining  filmwise  condensation  by  oxidizing  the 
tubes  with  a  special  solution  consisting  of  equal  amounts  of  sodium 
hydroxide  and  ethyl  alcohol  heated  to  about  80°  C.  The  oxidation  left  a 
very  thin  coating  on  the  tube  with  negligible  thermal  resistance. 
Although  Poole  had  troubles  in  maintaining  complete  filmwise  con¬ 
densation.  he  concluded  that  the  optimum  spacing  for  steam  conden¬ 
sation  on  a  horizontal  tube  with  rectangular  fins  is  around  1.5  mm. 
Poole  also  found  that  the  enhancement  of  the  finned  tube  over  the 
smooth  tube  was  greater  than  the  enhancement  in  surface  area  alone. 

Georgiadis  [Ref.  5],  in  September  1984,  improved  the  tube 
cleaning  and  oxidizing  procedures  and  verified  the  optimum  fin  spac¬ 
ing  for  steam  as  1.5  mm  for  rectangular  finned  tubes  with  various  fin 
thicknesses  and  heights.  He  found  that  the  heat- transfer  enhancement 
does  not  strongly  depend  on  fin  thickness,  though  a  thickness  of  0.75 
to  1.0  mm  thickness  was  slightly  better  than  the  average.  Although 
increased  fin  height  improved  enhancement,  increasing  fin  height 
above  1.0  mm  increased  the  heat-transfer  enhancement  less  than  the 
percentage  surface  area  gained.  Flook  [Ref.  9]  tested  tubes  with  trian¬ 
gular,  trapezoidal,  parabolic,  and  rectangular  shape  fins.  The  parabolic 
tubes  had  superior  heat-transfer  performance.  Flook  also  tested  for  the 
effect  of  tube  material  and  vapor  velocity.  He  found  that  an  increase  in 
vapor  velocity  from  2  to  8  m/s  increased  the  heat-transfer  rates  as 
much  as  ten  percent.  He  compared  the  results  of  a  stainless  steel  tube 
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with  a  copper  tube  and  confirmed  that  materials  with  high  thermal 
conductivities  will  enhance  the  heat-transfer  performance. 

In  March  1986.  Mitrou  (Ref.  10]  tested  wire-wrapped  tubes  and 
tubes  with  rectangular,  trapezoidal,  triangular,  and  parabolic  shape 
fins.  The  parabolic  fins  had  the  best  heat-transfer  performance.  Mitrou 
confirmed  the  effects  of  thermal  conductivity  discussed  by  Flook  by 
testing  tubes  made  out  of  copper,  aluminum,  copper-nickel,  and 
stainless  steel.  Cakan  [Ref.  11],  in  December  of  1986,  continued  test¬ 
ing  with  steam  and  attempted  to  improve  the  vapor-side  heat-transfer 
coefficient  by  attaching  a  porous  drainage  strip  to  the  underside  of  the 
finned  tube.  The  drainage  strip  was  designed  to  remove  the  thick 
condensate  film  on  the  lower  portion  of  the  tube  and  thus  enhance  the 
performance  of  the  lower  portion.  Cakan  found  that  the  use  of  drainage 
strips  significantly  enhanced  the  vapor-side  coefficient. 

In  June  1987,  Zebrowski  [Ref.  12]  was  the  first  student  to  test  a 
fluid  other  than  steam.  R-113  was  selected  because  of  its  low  surface 
tension.  Oxidizing  the  tube  was  not  required  to  maintain  complete 
filmwise  condensation  for  R-113.  Zebrowski's  research  resulted  in  an 
optimum  fin  spacing  between  0.25  and  0.5  mm.  He  showed  that  the 
Beatty-Katz  relationship  for  the  outside  heat-transfer  coefficient 
agreed  with  his  experimental  results  to  within  ±10  percent  for  inter- 
fin  spacings  greater  than  1.0  mm.  As  the  interfin  spacing  was 
decreased  from  1 .0  mm,  the  Beatty-Katz  relationship  overpredicted 
experimental  results.  Zebrowski  concluded  that  this  overprediction 
was  caused  by  condensate  retention  on  the  lower  portion  of  the  tube. 
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Zebrowski  also  measured  the  local  vapor-side  heat  transfer  coefficient 
for  two  finned  tubes.  This  was  accomplished  by  systematically  insulat¬ 
ing  the  upper  portions  of  the  tube  and  using  the  modified  Wilson  plot 
technique.  He  found  that  the  local  enhancement  at  the  top  of  the  tube 
was  approximately  twice  the  average  enhancement.  Lester  [Ref.  13],  in 
September  1987,  with  steam  as  the  working  fluid,  tested  for  the  local 
vapor-side  coefficient  around  a  finned  tube.  As  with  Zebrowski,  he 
concluded  that  at  the  top  of  the  tube,  the  local  coefficient  was  approx¬ 
imately  twice  the  average  value  around  the  tube. 

Table  1.1  shows  a  summary  of  the  NPS  research  program.  All  of 
the  previous  research  was  for  tubes  with  a  root  diameter  of  19.05  mm. 
Since  tube  diameter  may  play  a  significant  role  in  determining  the 
condensate  retention  angle  and  heat  and  fluid  flow  processes,  it  is 
important  that  this  aspect  of  the  problem  be  studied.  Also,  a  fluid  with 
intermediate  surface  tension  properties,  such  as  ethylene  glycol, 
should  be  tested  to  further  study  the  effect  of  fluid  properties  r  :  the 
condensation  process. 

E.  OBJECTIVES 

The  main  objectives  of  this  thesis  were: 

1.  to  measure  the  effect  of  root  diameter  on  the  enhancement  of 
finned  tubes,  and 

2.  to  modify  the  apparatus  to  enable  the  use  of  ethylene  glycol  as  a 
working  fluid. 
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Dec  86  Cakan  Steam  •  Found  enhancement  can  be  significantly  improved  by 

attaching  a  porous  drainage  strip  on  the  underside  of  a 
finned  tube. 
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n.  LITERATURE  SURVEY 


A.  INTRODUCTION 

Until  recently,  the  use  of  horizontal  low  integral-fin  tubes  for  heat 
transfer  with  high  surface  tension  working  fluids  (i.e.,  water)  was  con¬ 
sidered  impractical.  A  visual  inspection  of  the  of  the  horizontal  tube 
revealed  a  thick  film  of  retained  condensate  on  the  lower  portion  of 
the  tube.  Owing  to  this  thick  film,  the  heat- transfer  performance  in 
this  lower  portion  of  the  tube  was  expected  to  be  negligible.  However, 
Wanniarachchi.  et  al.  (Ref.  14]  and  Yau,  et  al.  [Ref.  4]  have  shown  con¬ 
siderable  enhancements  even  for  fully  flooded  tubes.  Therefore,  a 
complete  understanding  of  the  actual  physical  phenomena  that  take 
place  in  the  flooded  portion  is  essential  in  predicting  the  heat-transfer 
performance  of  finned  tubes  undergoing  filmwise  condensation. 

B.  CONDENSATE  RETENTION  ANGLE 

The  condensate  retention  angle  has  been  defined  as  the  angle 
measured  from  the  bottom  of  a  horizontal  finned  tube  to  a  point 
around  the  tube  circumference  where  the  condensate  film  between 
fins  abruptly  thins.  This  angle  is  graphically  depicted  in  Figure  2.1. 
Condensate  retention  angle  is  dependent  on  fluid  properties  and  fin 
geometry. 

Katz,  et  al.  [Ref.  15)  were  the  first  to  investigate  the  condensate 
retention  angle  and  to  develop  an  equation  for  predicting  its  behavior. 


10 


External  Diameter  of  fins 


Figure  2. 1  Schematic  of  Condensate  Retention  Angle 

on  Finned  Tubes 


Their  initial  tests  were  conducted  on  two  disks  separated  by  a  washer 
to  simulate  two  adjacent  fins.  Later,  tubes  with  fin  densities  of  276  to 
984  fins  per  meter  (fpm)  and  fin  heights  of  1.2  to  5.7  mm  were  tested. 
The  amount  of  condensate  retained  on  the  tube  and  the  retention 
angle  were  measured  by  weighing  and  by  visual  inspection  with  a 
cathetometer.  The  data  were  taken  in  a  static  condition  (i.e.,  no  con¬ 
densation  occurring)  with  water,  aniline,  acetone,  and  carbon  tetra¬ 
chloride.  The  surface  tension  for  each  of  these  fluids  was  measured  by 
the  pendant  drop  method  and  a  capillary  tube.  Katz  found  that  reten¬ 
tion  angle  was  dependent  on  fin  geometry  and  fluid  properties.  The 
ratio  of  the  surface  tension  to  the  difference  of  condensate  and  vapor 
densities  was  very  significant  and  explained  why  retention  angle  varied 
from  one  fluid  to  another  for  a  given  tube.  Since  the  vapor  density  is 
much  smaller  than  the  condensate  density,  it  was  ignored,  thereby 
making  the  ratio  of  surface  tension  to  condensate  density  the  major 
contributor  from  the  working  fluid.  Equation  2.1  shows  the  relation¬ 
ship  developed  by  Katz,  et  al.  (Ref.  15). 


V  a 

4Df  -  2  Dr  -  2  s 

(180) 

sin1!'  P,g 

f(tff-D2r)5 

V980 ) 

where 

'F  =  condensate  retention  angle  (degrees), 
o  =  surface  tension  (N/m), 
g  =  gravitational  constant  (9.806  m/s2), 
pi  =  condensate  density  (kg/m3). 
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Dr  =  root  diameter  (m), 

Df  =  fin  diameter  (m).  and 
s  =  interfin  spacing  for  rectangular  fin(m). 

This  equation  shows  that,  keeping  the  fin  height  and  spacing  constant, 
an  increase  in  root  diameter  would  result  in  a  decrease  in  condensate 
retention  angle. 

Until  an  investigation  by  Rudy  and  Webb  [Ref.  161  in  1981,  the 
retention  angle  had  been  determined  only  under  static  conditions. 
Rudy  and  Webb  measured  liquid  retention  angles  by  visual  sighting 
through  a  cathetometer.  They  measured  the  retention  angle  for  water, 
R-ll,  and  n-pentane  under  static  conditions  and  their  results  com¬ 
pared  favorably  with  the  model  developed  by  Katz,  et  al.  (Equation  2.1). 
The  retention  angles  for  R-ll  and  n-pentane  were  also  measured 
under  dynamic  conditions.  Integral-fin  tubes  with  748,  1024,  1378 
fpm,  a  “Thermoexcel-CTtube  with  1429  fpm,  and  a  pin-fin  tube  with 
1378  fpm  were  examined  under  static  and  dynamic  conditions.  They 
confirmed  the  Katz,  et  al.  [Ref.  151  conclusion  that  retention  angle 
increased  with  the  ratio  of  surface  tension  to  condensate  density,  and 
concluded  that  the  condensate  retention  angle  was  not  significantly 
different  under  dynamic  and  static  conditions. 

In  1982,  Rifert  [Ref.  17)  developed  an  equation  to  predict  con¬ 
densate  retention  angle  by  modeling  the  rise  height  between  fins  to 
capillary  rise  height  along  vertical  plates.  Again,  retention  angle  was 
strongly  influenced  by  the  ratio  of  surface  tension  to  condensate 
density.  Equation  2.2,  developed  by  Rifert,  demonstrates  that  an 
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increase  in  root  diameter  will  decrease  the  condensate  retention 
angle 

'F  =  cos-1 

where 

Ap  =  profile  area  of  fin  (m2), 

Pw  =  wetted  perimeter  (m),  and 
p  =  fin  pitch  (m). 

Rudy  and  Webb  [Ref.  181  developed  a  general  equation  to  predict 
retention  angle  for  finned  tubes  of  various  fin  geometries.  Their 
working  fluids  of  water,  R-ll,  n-pentane,  and  acetone  were  used  with 
finned  tubes  having  fin  densities  of  748  to  1417  fpm.  Their  static  test 
consisted  of  a  finned  tube  being  placed  in  a  shallow  pool  of  working 
fluid.  The  rise  height  of  the  working  fluid  in  the  interfin  spaces  was 
compared  to  an  unrolled  tube  of  the  same  fin  geometry.  The  unrolled 
tube  was  fabricated  by  slicing  a  tube  section  and  carefully  peeling  the 
tube  open  into  a  flat  plate.  The  rise  heights  were  the  same.  When  the 
finned  tube  and  the  unrolled  tube  were  removed  from  the  shallow 
pool,  the  rise  heights  remained  unchanged.  From  these  observations, 
Rudy  and  Webb  deduced  that  the  retention  angle  was  caused  by  a  bal¬ 
ance  of  surface  tension  forces  and  the  weight  of  condensate.  A  simple 
force  balance  on  the  retained  condensate  led  to  the  following  equation: 


1  - 


2o(pw~P) 

P,gDfAp 


(2.2) 


T  =  cos 


2  o(2  e  - 1) 
Dfp,ges 


(2.3) 
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where 


e  =  fin  height  (m),  and 
t  =  fin  thickness  for  rectangular  fin  (m). 

Rudy  and  Webb  found  that  Equation  2.3  predicted  condensate  reten¬ 
tion  angle  to  within  ±10  percent.  Again,  increasing  the  root  diameter 
would  lower  the  condensate  retention  angle. 

Honda,  et  al.  [Ref.  19]  did  a  photographic  study  on  dynamic  and 
static  condensate  retention  angles  for  methanol  and  R-113.  They  con¬ 
firmed  the  conclusion  of  Rudy  and  Webb  that  static  and  dynamic 
retention  angles  were  approximately  the  same.  They  experimented  on 
finned  tubes  with  and  without  porous  drainage  strips.  Their  measure¬ 
ments  with  no  strips  were  in  agreement  with  previous  experimental 
data,  and  they  developed  the  following  equation: 


¥  =  cos-1 


4 ocos p 
DfP,gs 


(2.4) 


where 

P  =  fin-tip  half  angle. 

In  1985,  Rudy  and  Webb  [Ref.  20}  extended  their  prediction  of 
condensate  retention  angle  for  various  fin  geometries  by  modifying 
Equation  2.3.  They  tested  finned  tubes  with  fin  densities  of  748  to 
1378  fpm,  one  spine  fin  tube  with  1378  fpm,  and  a  "Thermoexcel-CT 
tube.  They  used  R-ll,  water,  and  n-pentane  as  the  working  fluids. 
Again,  Rudy  and  Webb  concluded  that  the  retention  angle  increases 
with  the  ratio  of  surface  tension  to  condensate  density  and  fin  height 
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but  decreases  as  tube  diameter  increases.  They  developed  the  follow¬ 
ing  expression  for  a  trapezoidal  shape  fin  shown  in  Figure  2.2. 


¥  =  cos'1 1  - 


2q(tt  +  2e-tb) 

DfP1g(tte"sbe-etb)J  * 


(2.5) 


where 

tt  =  fin  thickness  at  tip  (m). 
tb  =  fin  thickness  at  base  (m),  and 
Sb  =  interfin  spacing  at  base  (m). 

For  a  rectangular-shape  fin.  Equation  2.5  reduces  to  the  following: 


=  cos 


-i ,  4o 

DfPigs 


(2.6) 


This  equation  is  identical  to  Equation  2.4  of  Honda  when  (5=0 
degrees. 


Figure  2.2  Trapezoidal  Fin  Used  by  Rudy  and  Webb  [Ref.  20] 
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Once  again.  Equations  2.5  and  2.6  predict  that  condensate  reten¬ 
tion  angle  will  decrease  with  increasing  tube  diameter.  Rudy  and  Webb 
completed  static  retention  angle  tests  on  19  mm  and  25  mm  diameter 
tubes  with  fin  densities  of  1024  fpm.  They  observed  that  the  rise 
heights  were  equal  and  therefore  the  larger-diameter  tube  retained 
less  condensate  and  had  a  smaller  retention  angle. 

Masuda  and  Rose  (Ref.  21]  found  that  condensate  was  retained  not 
only  in  the  flooded  portion  of  the  tube.  Their  photographic  study 
revealed  that  condensate  was  also  retained  at  the  base  of  the  fin  in  the 
upper  portion  of  the  tube  usually  referred  to  as  the  unflooded  portion. 
The  film  thickness  around  the  base  decreased  with  circumferential 
position  around  the  tube.  They  therefore  postulated  that  more  surface 
area  of  the  finned  tube  was  insulated  by  the  condensate  than  previously 
expected.  Masuda  and  Rose  (Ref.  211  concentrated  their  research  on 
the  profile  of  the  condensate  in  the  interfin  spacing.  They  defined  four 
separate  retention  conditions  based  on  the  meniscus  profile  in  the 
interfin  spaces  and  along  the  sides  of  the  fin.  These  retention  condi¬ 
tions  are  listed  below  and  shown  pictorially  in  Figure  2.3. 

1.  “...the  interfin  space  is  just  filled  by  the  meniscus  but  the  fin 
flanks  are  not  wholly  wetted."  [Ref.  21]  (Figure  2.3  b(2)) 

2.  “...where  the  whole  of  the  flank  is  just  wetted  and  for  which  the 
liquid  film  at  the  center  of  the  interfin  spacing  has  finite  thick¬ 
ness."  (Ref.  21)  (Figure  2.3  b(3),  defined  as  fully  flooded) 

3.  “...when  the  fin  flank  is  just  completely  wetted  but  the  interfin 
space  is  not  wholly  wetted."  [Ref.  21]  (Figure  2.3  c(2)) 

4.  “...the  whole  of  the  interfin  space  is  Just  wetted  and  for  which  the 
contact  angle  at  the  fin  tip  is  non-zero."  [Ref.  21]  (Figure  2.3  c(3)) 
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Masuda  and  Rose  developed  separate  expressions  for  each  listed 
flooding  condition.  Condition  three  was  considered  the  fully  flooded 
condition  and,  in  this  situation,  their  expression  reduces  to  Equation 
2.6  for  rectangular  fins. 

Honda,  et  al.  [Ref.  22],  in  1987,  developed  two  equations  for  the 
prediction  of  condensate  retention  angle.  One  expression  for  retention 
angle  was  used  when  the  fin  spacing  exceeded  twice  the  fin  height 
and  the  other  when  the  fin  spacing  was  less  than  twice  the  fin  height. 
The  latter  case,  for  rectangular  finned  tubes,  reduces  to  Equation  2.6. 


C.  HEAT-TRANSFER  PERFORMANCE  ON  FINNED  TUBES 

In  1948,  Beatty  and  Katz  [Ref.  23)  were  the  first  to  develop  a 
model  to  predict  the  outside  condensing  heat-transfer  coefficient  for  a 
horizontal  finned  tube.  They  divided  the  finned  tube  into  two  geomet¬ 
ric  regions.  The  fin  flank  was  considered  to  be  a  vertical  plate  and  the 
interfin  spacing  to  be  a  smooth  horizontal  tube.  By  combining  the 
equations  of  Nusselt  for  vertical  plates  and  smooth  horizontal  tubes, 
Beatty  and  Katz  expressed  the  outside  heat-transfer  coefficient  as 
shown  in  Equation  2.7. 


hBK  =0.689F 


029 


A  ^  no.25 ) +  f -3tj  f- 
** o  n 


o  L 


025 


(2.7) 


where 

k®Plghfe 

F=|_  H.ATv, 


(2.8) 
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and 


(2.9) 


hsK  =  Beatty-Katz  outside  heat-transfer  coefficient(W/m2  •  K), 
hfg  =  specific  enthalpy  of  vaporization  (J/kg), 
ki  =  condensate  thermal  conductivity  (W/m  •  K). 
pi  =  condensate  dynamic  viscosity  (kg/m*s), 

T]  f  =  fin  efficiency. 

ATVS  =  temperature  drop  across  the  condensate  film  (K). 

Ar  =  surface  area  of  smooth  tube  (m2), 

Af  =  surface  area  of  fin  (m2),  and 

Aq  =  total  external  surface  area,  Ar  +  nf  Af,  (m2). 

The  coefficient  of  0.689  is  five  percent  less  than  the  expected  value  of 
0.728  from  the  Nusselt  equation  for  a  smooth  horizontal  tube.  This 
coefficient  resulted  from  tests  carried  out  on  six  single  horizontal 
finned  tubes.  These  tubes  had  fin  densities  of  up  to  630  fpm.  Their 
working  fluids  included  methyl  chloride,  sulfur  dioxide,  refrigerant- 
22,  propane,  n-butane,  and  n-pentane.  Surface  tension  for  these  fluids 
varied  from  0.005  to  0.025  N/m.  They  found  that  the  computed  results 
from  Equation  2.7  and  experimental  data  agreed  within  ten  percent. 
Notice  that  Equation  2.7  shows  that  the  heat-transfer  coefficient 
decreases  with  increasing  tube  diameter. 

Although  the  use  of  Nusselt’s  equations  was  ingenious,  those 
equations  forced  Beatty  and  Katz  to  apply  Nusselt’s  assumptions  to 
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their  model.  The  most  important  assumption  was  that  gravity  forces 
alone  controlled  condensate  drainage.  No  surface  tension  forces 
existed  and  no  condensate  was  retained  between  fins.  Because  Beatty 
and  Katz  [Ref.  23]  only  tested  tubes  with  low  fin  densities  and  used 
relatively  low  surface  tension  working  fluids,  their  equation  seemed 
promising.  Rudy  and  Webb  [Ref.  16]  found  that  the  Beatty  and  Katz 
model  overpredicted  the  results  for  fluids  with  a  higher  surface  ten¬ 
sion  to  condensate  density  ratio.  They  attempted  to  use  the  Beatty  and 
Katz  relationship  for  only  the  unflooded  region,  which  led  to  Equation 
2.10. 


where  Iibk  is  given  by  Equation  2.7  and 

h*BK  =  modified  Beatty-Katz  condensing  coefficient  (W/m2  •  K). 

This  was  the  first  research  to  include  the  effect  of  condensate  reten¬ 
tion  in  heat  transfer.  They  reasoned  that  surface  tension  forces  caused 
the  retention  of  the  condensate  on  the  lower  portion  of  the  tube  to  act 
as  an  insulator.  This  equation  resulted  in  10-  to  50-percent  underpre¬ 
dictions  for  R-ll,  n-pentane,  and  water. 

Owen,  et  al.  [Ref.  24],  in  1983,  concluded  that  the  margin  of  error 
in  the  Rudy  and  Webb  Equation  resulted  from  the  assumption  of  no 
heat  transfer  in  the  flooded  region.  Therefore,  to  improve  this  model, 
Owen,  et  al.  allowed  for  one-dimensional  heat  conduction  through  the 
fins  and  condensate  film  in  the  flooded  region.  The  condensate  thick¬ 
ness  was  presumed  to  be  equal  to  the  fin  height  in  the  flooded  region. 
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This  model,  however,  showed  only  a  slight  improvement  over  the 
Rudy  and  Webb  model. 

It  became  apparent  that  condensation  on  a  finned  tube  was  a  very 
complex  phenomenon  with  a  number  of  controlling  parameters  not 
considered  in  previous  models.  Variables  not  considered  were  three- 
dimensional  condensate  flow,  surface  tension  forces,  wall  conduction 
effects,  condensate  film  thickness  variations,  and  vapor  velocities  [Ref. 
25[.  Previous  models  considered  only  the  insulating  effect  of  the  con¬ 
densate  film  in  the  flooded  region.  Gregorig  [Ref.  26),  as  early  as  in 
1954,  proved  that  surface  tension  also  produces  pressure  gradients 
that  thin  the  condensate  film  around  the  fin  tips  and  along  the  fin 
sides. 

In  1987,  Webb,  et  al.  [Ref.  27]  re-examined  the  Beatty-Katz  model 
and  pointed  out  that  on  finned  tubes  surface  tension,  not  gravity,  was 
the  prime  mover  of  condensate.  Therefore,  the  Nusselt  analysis  was 
not  valid  for  finned  tubes.  They  therefore  divided  the  finned  tube  into 
two  regions:  the  flooded  and  the  unflooded  portions.  Equation  2.11 
was  used  to  compute  the  average  outside  heat-transfer  coefficient 


hq  = 


180-^" 

180 


Ar  u 
hY  +  hf1 


(2.11) 


where 

h  =  condensing  coefficient  based  on  A^f  (W/m2  •  K), 
hf  =  condensing  coefficient  of  the  fin  (W/m2  •  K), 
hh  =  condensing  coefficient  for  smooth  horizontal  tube  (W/m2  •  K), 
hb  =  condensing  coefficient  for  flooded  portion  (W/m2  •  K), 
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Adf  =7tDf  L  (m2). 

L  =  length  of  tube  (m),  and 

rj  =  surface  efficiency,  (1  -  (1  -  rjf)  Ar/(Ar+Af)). 

The  condensing  coefficient  for  the  fin  (hf)  was  determined  by 
assuming  that  surface  tension  forces  induced  pressure  gradients  along 
the  condensate  film  surface.  They  assumed  a  linear  pressure  gradient 
from  fin  tip  to  fin  root.  A  heat  balance  on  the  condensate  yielded  the 
average  heat-transfer  coefficient  for  the  fin.  To  compute  the  heat- 
transfer  coefficient  for  the  horizontal  smooth  tube  portion  (hh)  in  the 
interfin  spaces,  an  iterative  process  which  would  take  into  account  the 
additional  condensate  from  the  fin  flanks  was  used.  The  condensing 
coefficient  in  the  flooded  region  (hb)  was  derived  by  use  of  two- 
dimensional  conduction  of  the  fin-condensate  composite  structure. 
They  found  that  in  the  flooded  region,  the  heat-transfer  rate  only 
accounted  for  0.2  percent  of  the  total  rate  for  R-ll.  For  steam  with  a 
retention  angle  of  54  degrees,  the  heat  transfer  rate  in  the  flooded 
region  was  only  1.6  percent  of  the  total  rate.  These  results  suggested 
that  the  flooded  portion  of  the  tube  could  be  practically  ignored. 
Notice  that  these  calculated  results  are  in  disagreement  with  the 
experimentally  observed  results  mentioned  earlier  that  even  fully 
flooded  tubes  gave  considerable  enhancements.  Therefore,  it  is  possi¬ 
ble  that  actual  processes  with  fully  flooded  tubes  are  far  more  complex 
than  can  be  explained  by  one-  or  two-dimensional  heat  conduction 
alone.  By  ignoring  the  heat  transfer  in  the  flooded  portion,  Webb,  et  al. 
(Ref.  27]  pointed  out  the  importance  of  the  retention  angle.  By 
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decreasing  the  retention  angle,  the  heat  transfer  would  significantly 
increase  as  shown  by  Equation  2.11.  With  this  model,  Webb,  et  al.  [Ref. 
27]  found  that  experimental  and  theoretical  results  for  R-ll  agreed  to 
within  ±25  percent.  Disappointed  with  these  results,  Webb,  et  al.  [Ref. 
27]  used  the  theory  of  Adamek  [Ref.  28]  to  predict  a  new  value  for  the 
condensing  coefficient  on  the  unflooded  fin  flanks.  This  latter  method 
provided  results  within  ±20  percent  of  their  experimental  data  for 
R-ll. 

In  1987,  Honda,  et  al.  [Ref.  29]  introduced  a  prediction  model 
that  accounted  for  surface-tension-induced  pressure  gradients,  con¬ 
densate  retention  angle,  heat  transfer  through  the  flooded  region,  and 
fin  geometry.  The  tube  was  divided  into  flooded  and  unflooded  regions 
and  the  average  Nusselt  number  and  wall  temperature  were  computed 
for  each  region. 


Nud  = 


, Tl  u(l  -Tw)(l  -  fr)  +  Nu  df T1  f(l  - T 

(l-T^a-'n  +  U-Tj* 


(2.12) 


where  the  dimensionless  temperature  (T)  and  retention  angle  (^J 
were  expressed  as: 


(2.13) 


and 


(2.14) 
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Nud  =  average  Nusselt  number  for  finned  tube, 

Nudu  =  average  Nusselt  number  for  unflooded  portion, 

Nudf  =  average  Nusselt  number  for  flooded  portion, 

Tc  =  average  coolant  temperature, 

Ts  =  saturation  temperature, 

T 

=  dimensionless  average  wall  temperature  for  flooded  por¬ 
tion.  and 

?wu  =  dimensionless  average  wall  temperature  for  unflooded 
portion. 

Equation  2.6  was  used  to  solve  for  the  condensate  retention  angle.  In 
Equation  2.12,  the  retention  angle  plays  a  significant  role  in  deter¬ 
mining  the  heat  transfer.  The  average  Nusselt  numbers  were  used  in 
Equation  2.12  because  Honda,  et  al.  [Ref.  29]  allowed  the  surface- 
tension-induced  pressure  gradients  to  change  around  the  circumfer¬ 
ence  of  the  tube.  A  ±20  percent  error  was  found  between  this  model 
and  experimental  results  for  22  tubes  and  11  fluids.  Marto,  et  al.  [Ref. 
30]  compared  this  model  with  their  R-113  experimental  results  and 
found  the  model  to  be  generally  conservative  with  a  ±20  percent  error. 
The  model  did  underpredict  for  very  small  fin  spacings. 

D.  SUMMARY 

From  these  previous  findings,  it  is  apparent  that  the  condensate 
retention  angle  can  be  easily  predicted.  Equation  2.6  is  the  generally 
accepted  expression  for  condensate  retention  angle.  From  this  equa¬ 
tion,  it  is  easily  seen  that  an  increase  in  the  tube  root  diameter  will 
decrease  the  condensate  retention  angle.  Although  a  great  deal  of 
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effort  has  been  applied  to  develop  an  accurate  heat-transfer  prediction 
model,  one  does  not  exist.  Marto.  et  al.  (Ref.  30]  found  that  the  Honda, 
et  al.  (Ref.  29]  model  appears  to  be  the  best,  showing  an  agreement  to 
within  ±20  percent  of  the  experimental  data.  This  percent  error  was 
far  better  over  the  entire  range  of  fin  spacings  tested  than  obtained 
with  the  Webb,  et  al.  (Ref.  27]  and  the  Beatty-Katz  (Ref.  23]  models. 
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m.  APPARATUS  AND  TUBES  TESTED 


A.  DESCRIPTION  OF  APPARATUS 

The  test  apparatus,  as  shown  in  Figure  3.1,  consisted  of  a  boiler,  a 
test  section,  glass  piping,  and  an  auxiliary  condenser  with  a  purging 
system  to  remove  non-condensing  gases.  The  boiler  was  fabricated 
from  a  0.3048  m  (12  inch)  diameter  Pyrex  glass  section  with  ten 
4-kW,  440-V,  Watlow  immersion  heaters.  Vapor  flowed  upward  from 
the  boiler  into  a  2.44  m-long  section  of  152  mm-diameter  glass  piping 
before  making  a  180-degree  turn  toward  the  test  section  1.52  m 
below.  The  test  section  (see  Figure  3.2)  was  fabricated  from  stainless 
steel  with  nylon  and  Teflon  fittings  to  support  the  test  tube  and  to 
connect  cooling  water.  Cooling  water  to  the  test  tube  was  pumped 
from  a  water  sump  and  was  varied  by  a  throttle  valve  at  the  inlet  of  a 
flowmeter.  A  continuous  flow  of  tap  water  was  supplied  to  the  water 
sump.  A  thermocouple,  placed  just  above  the  test  tube,  and  a  mano¬ 
meter  provided  accurate  vapor  temperature  and  pressure  readings, 
respectively.  A  view  port  in  the  test  section  allowed  the  condensation 
process  to  be  observed  and  photographed.  The  auxiliary  condenser 
consisted  of  a  copper  coil  within  Pyrex  glass  piping.  Cooling  water 
supply  was  tap  water  that  was  controlled  by  a  throttle  valve  at  the  inlet 
of  a  flowmeter.  Modifications  made  to  the  auxiliary  condenser  are 
discussed  later  in  this  chapter. 
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Figure  3.2  Schematic  of  Test  Section  (Insert  Removed) 
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The  purging  system  consisted  of  copper  tubing,  a  small  copper 
coil,  a  Plexiglas  cylinder,  and  a  vacuum  pump.  Gases  drawn  by  the 
pump  left  the  secondary  condenser  via  the  copper  tubing  and  flowed 
into  the  copper  coil.  This  coil  was  located  in  the  water  sump  that  pro¬ 
vided  cooling  water  to  the  test  tube.  The  heat  exchange  in  the  coil 
condensed  any  residual  working  fluid  vapor.  The  non-condensing 
gases  and  residual  working  fluid  were  collected  in  a  small  cylinder  in 
the  water  sump.  The  liquid  collected  at  the  bottom  of  the  cylinder  was 
drained  at  the  end  of  the  data  run.  The  non-condensing  gases  (usually 
air)  were  suctioned  to  the  top  of  the  container  and  drawn  into  the 
pump  suction  (Figure  3.3).  A  more  detailed  description  of  the  appara¬ 
tus  has  been  provided  in  the  references  [Refs.  5,  8,  9,  11,  13). 

B.  MODIFICATIONS  TO  APPARATUS 

Modifications  were  made  to  the  original  apparatus  to  improve  sys¬ 
tem  integrity  and  to  meet  the  objectives  of  this  thesis.  Prior  to  these 
modifications,  cooling  water  to  the  secondary  condenser  was  con¬ 
tained  in  two  helically  wound  coils  made  of  3/8-inch  copper  tubing.  A 
4-inch  coil  was  located  inside  a  5  1 /2-inch  coil.  Each  coil  was  18 
inches  in  length.  Preliminary  tests  with  ethylene  glycol  at  a  saturation 
temperature  of  128°  C  indicated  that  the  water  flow  rate  through 
these  coils  was  not  large  enough.  When  proper  water  flow  rate  was 
provided  to  achieve  a  desired  system  pressure,  slugs  of  liquid  and 
steam  were  observed  to  be  flowing  out  of  the  cooling  coils,  indicating 
that  boiling  of  the  water  was  occurring.  The  steam  produced  created 
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vapor-blocking  inside  the  tubes,  thus  almost  stopping  the  water  flow 
intermittently.  When  this  occurred,  the  system  pressure  started  to 
increase.  Once  the  tube  was  entirely  blanketed  with  vapor,  water 
rushed  again,  providing  cooling.  This  cycle  repeated  while  the  system 
pressure  experienced  uncontrollable  fluctuations.  To  avoid  this  prob¬ 
lem,  a  new  5-inch  diameter  helical  coil  was  therefore  fabricated  from 
5/8-inch  copper  tubing.  The  new  coil  was  18  inches  in  length  and 
tested  satisfactorily  for  all  three  working  fluids.  The  larger  inner  dia¬ 
meter  of  this  tube  was  adequate  in  preventing  any  vapor  blocking  by 
allowing  a  larger  cooling  water  mass  flow  rate. 

Since  the  new  configuration  called  for  a  single  coil  over  two,  a 
new  base  plate  for  the  secondary  condenser  was  needed.  The  fitting 
for  the  gravity  drain  condensate  return  line  was  1/2  inch  vice  the  1- 
inch  fitting  out  of  the  dual  coil  base  plate.  A  single  1 /2-inch  stainless 
steel  tube  with  a  90  degree  bend  was  used  for  the  drain  line.  Previ¬ 
ously,  the  old  drain  line  reduced  from  1  inch  to  1/2  inch  by  using  two 
reducing  fittings.  This  drain  line  was  a  constant  source  of  vacuum 
leaks  and  its  location  in  the  system  made  repairs  very  difficult.  The 
installation  of  the  single  1  /2-inch  line  improved  system  integrity  by 
eliminating  a  possible  source  of  vacuum  leaks. 

C.  TUBES  TESTED 

One  smooth  tube,  18  tubes  having  rectangular  section  fins,  and  a 
special  tube  identical  to  the  one  used  by  Masuda,  et  al.  [Ref.  3]  of 
Queen  Mary  College,  University  of  London,  were  to  be  tested  with 
steam,  R-113,  and  ethylene  glycol  as  the  working  fluids.  All  tubes  were 
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made  of  copper  and  the  18  finned  tubes  consisted  of  three  families 
having  different  tube  root  diameters.  The  tube  diameter  was  the  char¬ 
acteristic  common  to  each  family.  The  tubes  with  the  12.7  mm  diame¬ 
ter  (also  referred  to  as  root  diameter)  were  considered  the  small 
tubes,  while  the  19.05  mm  and  the  25  mm  were  medium  and  large 
tubes,  respectively.  Each  tube  matched  one  other  tube  in  a  different 
family  with  all  geometric  parameters,  with  the  exception  of  the 
diameter.  All  tubes  had  a  fin  height  of  1 .0  mm  and  a  fin  thickness  of 
1.0  mm.  The  comparison  tube  will  be  referred  to  as  the  QMC  tube. 
Tube  parameters  are  listed  in  Table  3.1.  Figure  3.4  shows  a  photo¬ 
graph  of  the  small-diameter  tube  family.  Figure  3.5  shows  a  photo¬ 
graph  of  tubes  with  a  1.5  mm  fin  spacing,  while  Figure  3.6  shows  a 
photograph  of  the  inserts  used  during  this  investigation. 

Various  problems  were  encountered  while  testing  and  will  be  dis¬ 
cussed  in  Chapter  V.  Only  preliminary  data  were  taken  for  the  small 
and  QMC  tube  using  ethylene  glycol.  After  completing  R-113  and 
ethylene  glycol  data,  tube  10  was  mistakenly  modified  by  blasting  its 
surface  with  glass  beads  in  an  attempt  to  remove  the  oxide  film.  Close 
inspection  of  the  tube  surface  area  revealed  that  the  local  geometry 
around  the  tips  of  the  fins  had  been  significantly  altered.  Therefore,  a 
tube  of  the  same  dimensions  was  produced  to  replace  tube  10.  How¬ 
ever,  the  new  tube  was  manufactured  at  the  NPS  machine  shop  and  its 
copper  surface  was  slightly  different  than  the  original  tube  that  was 
commercially  manufactured. 
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TABLE  3.1 


DIMENSIONS  OP  TUBES  TESTED 


Tube 

Dr  (mm) 

Di  (mm) 

s  (mm) 

t  (mm) 

e  (mm) 

Smooth  tube 

1 

19.05a 

12.7 

- 

- 

- 

Small  tubes 

2 

12.7 

9.53 

0.25 

1.0 

1.0 

3 

12.7 

9.53 

0.5 

1.0 

1.0 

4 

12.7 

9.53 

1.0 

1.0 

1.0 

5 

12.7 

9.53 

1.5 

1.0 

1.0 

6 

12.7 

9.53 

2.0 

1.0 

1.0 

7 

12.7 

9.53 

4.0 

1.0 

1.0 

Medium  tubes 

8 

19.05 

12.7 

0.25 

1.0 

1.0 

9 

19.05 

12.7 

0.5 

1.0 

1.0 

10 

19.05b 

12.7 

1.0 

1.0 

1.0 

11 

19.05 

12.7 

1.5 

1.0 

1.0 

12 

19.05 

12.7 

2.0 

1.0 

1.0 

13 

19.05 

12.7 

4.0 

1.0 

1.0 

Large  tubes 

14 

25.0 

12.7 

0.25 

1.0 

1.0 

15 

25.0 

12.7 

0.5 

1.0 

1.0 

16 

25.0 

12.7 

1.0 

1.0 

1.0 

17 

25.0 

12.7 

1.5 

1.0 

1.0 

18 

25.0 

12.7 

2.0 

1.0 

1.0 

19 

25.0 

12.7 

4.0 

1.0 

1.0 

QMC  tube 

20 

12.7 

9.53 

1 

0.5 

1.59 

a  Dr  =  Do  for  smooth  tube 
b  Glass  beaded  tube 
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Figure  3.4  Photograph  of  Small-Diameter  Tube  Family 


Figure  3.5  Photograph  of  Tubes  With  1.5  mm  Spacing 


Figure  3.6  Photograph  of  Inserts 


D.  SYSTEM  INTEGRITY 

After  completing  the  modifications  to  the  apparatus,  it  was 
important  to  verify  system  integrity  prior  to  taking  data.  System  integ¬ 
rity  is  a  term  used  to  describe  the  air  tightness  of  the  closed-loop 
apparatus.  For  approximately  three  months,  the  system  had  remained 
idle  prior  to  this  research.  For  this  reason,  all  the  rubber  gaskets  were 
replaced.  A  vacuum  test  was  conducted  by  using  the  vacuum  system  to 
bring  the  pressure  in  the  apparatus  down  to  approximately  20  mm  of 
mercury.  The  apparatus  was  secured  and  the  vacuum  pump  turned  off. 
If  the  manometer  indicated  that  the  apparatus  had  maintained  a  vac¬ 
uum  within  2  mm  in  a  24-hour  period,  system  integrity  was  consid¬ 
ered  good.  This  test  was  repeated  between  data  runs  at  least  once  a 
week. 

Initially,  many  leaks  in  the  apparatus  were  found  around  the  fit¬ 
tings  of  the  condensate  drain  lines.  Three  valves  and  a  second  gasket 
were  replaced.  Some  leaks  were  located  by  placing  the  system  under 
vacuum  to  draw  in  Freon  gases.  Once  the  Freon  gases  were  in  the 
apparatus,  nitrogen  gas  was  used  to  slightly  pressurize  the  apparatus.  A 
Freon  detector  was  then  slowly  moved  around  the  outside  of  the  appa¬ 
ratus  to  detect  leaks.  This  method  was  successful  in  finding  the  last 
few  leaks. 

The  presence  of  non-condensing  gases  during  data  taking  was  a 
good  indication  of  a  vacuum  leak.  The  data  reduction  program  would 
ask  the  operator  if  non-condensing  gas  concentration  needed  to  be 
checked.  The  check  compared  the  actual  vapor  temperature  measured 
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by  the  thermocouple  in  the  vapor  space  upstream  of  the  test  tube  to 
the  calculated  saturation  temperature  for  the  apparatus  pressure.  Non¬ 
condensing  gases  were  checked  prior  to  the  first  data  set  accepted 
and  the  last  data  set  accepted.  If  the  check  prior  to  the  last  data  set 
showed  the  presence  of  non-condensing  gases,  the  entire  data  run  was 
discarded.  A  vacuum  test  was  conducted  after  non-condensing  gases 
were  detected  or  if  the  data  were  suspect.  Most  vacuum  leaks  were 
caused  by  a  bad  gasket  or  valve  on  the  secondary  condenser  or  by  the 
vacuum  system  being  inadvertently  left  open. 
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IV.  EXPERIMENTAL  AND  DATA  REDUCTION  PROCEDURES 

A.  EXPERIMENTAL  PROCEDURES 

Prior  to  testing,  a  soft  brush  and  tap  water  were  used  to  clean  the 
interior  and  exterior  of  each  tube.  Distilled  water  was  used  for  the  final 
rinse.  As  proven  by  Zebrowski  [Ref.  12],  while  testing  R-113,  an  oxi¬ 
dizing  coating  was  not  needed  to  assure  filmwise  condensation.  A 
preliminary  test  with  ethylene  glycol  indicated  that  there  was  no  dif¬ 
ference  between  the  results  with  and  without  a  coating.  For  this  rea¬ 
son,  it  was  determined  that  the  oxidizing  coating  was  not  required  for 
filmwise  condensation  with  ethylene  glycol. 

When  the  working  fluid  was  steam,  the  tube  had  to  be  coated  with 
a  thin  oxidation  layer  with  negligible  thermal  resistance  to  ensure 
filmwise  condensation.  Equal  volumes  of  sodium  hydroxide  and  ethyl 
alcohol  were  mixed  and  heated  to  about  80°  C  to  make  an  oxidizing 
solution.  This  solution  was  applied  while  the  tube  was  suspended  above 
a  pool  of  boiling  water.  To  establish  a  coating,  this  oxidizing  solution 
was  applied  every  ten  minutes  for  one  hour.  When  retesting  a  tube,  the 
solution  was  applied  every  five  minutes  for  20  minutes.  Once  a  tube 
was  coated,  it  was  rinsed  with  distilled  water.  Also,  with  steam,  an 
insert  was  placed  in  the  tube  to  enhance  the  internal  heat-transfer 
coefficient.  The  reason  for  this  procedure  is  explained  later  in  this 
chapter.  The  insert  was  placed  in  the  same  position  for  every  tube  (see 
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Figure  3.6)  Since  the  large-  and  medium-diameter  tubes  had  the  same 
inside  diameter,  the  same  insert  was  placed  in  them. 

When  the  tube  was  ready  for  testing,  it  was  installed  in  the  test 
section.  Installation  took  approximately  five  minutes.  The  tap  water 
supply  to  the  water  sump  was  opened  and  the  pumps  were  started.  An 
initial  flow  rate  of  20  percent  was  set  and  a  visual  inspection  for  leaks 
at  the  fittings  was  conducted.  The  boiler  control  panel  was  energized 
and  the  system  was  brought  up  to  operating  temperature  by  adjusting 
the  boiler  input  power,  adjusting  cooling  water  flow  to  the  secondary 
condenser,  and  operating  the  vacuum  system.  The  procedures  for 
placing  this  system  in  operation  were  outlined  by  Poole  [Ref.  8]. 

The  removal  of  non-condensing  gases  was  vital  to  ensure  reliable 
and  consistent  results.  The  Gibbs-Dalton  ideal  gas  mixture  relations, 
together  with  the  measured  vapor  temperature  and  the  temperature 
corresponding  to  the  measured  pressure,  were  used  to  compute  the 
non-condensing  gas  (believed  to  be  air)  concentration.  The  computed 
non-condensing  gas  molar  concentrations  were  found  to  be  ±0.5  per¬ 
cent.  This  reveals  a  non-condensing  gas  concentration  which  is  zero  to 
within  the  accuracy  of  the  measurements.  Once  a  tube  had  been 
installed  in  the  apparatus,  the  purge  system  was  operated  successfully 
to  remove  all  non-condensing  gases  before  storing  any  data.  When 
testing  with  steam,  the  vacuum  system  was  in  operation  with  the  inlet 
valve  opened  fully  for  low  pressure  and  opened  slightly  for  atmo¬ 
spheric  pressure.  Georgiadis  [Ref.  5]  found  that  operating  the  vacuum 
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pump  did  not  alter  experimental  results  but  maintained  filmwise  con¬ 
densation  conditions  longer. 

The  output  of  the  vapor  thermocouple,  which  was  continuously 
displayed  on  the  voltmeter,  was  used  in  obtaining  and  maintaining  the 
system  temperature  (hence  pressure).  Once  the  operating  tempera¬ 
ture  was  established,  sample  data  were  taken  to  check  for  non-con¬ 
densing  gases  and  desired  vapor  velocity.  For  example,  if  the  vapor 
velocity  was  below  the  desired  value,  the  boiler  power  and  the  flow  to 
the  secondary  condenser  cooling  water  were  increased.  Table  4.1 
shows  a  list  of  operating  conditions  for  each  fluid. 

TABLE  4.1 

OPERATING  CONDITIONS 


Fluid 

Thermocouple 

(microvolts) 

Temperature 

(°C) 

Pressure 
(mm  Hg) 

R-113 

1977 

48 

765 

Steam 

1977 

48 

85 

Steam 

4280 

100 

765 

Ethylene  Glycol 

5400 

128 

62 

For  each  data  run,  the  steady-state  condition  was  maintained  for 
about  30  minutes  before  any  taking  of  data.  At  each  data  point,  the 
test-tube  cooling  water  flow  rate  was  manually  entered.  When  R-113 
was  the  working  fluid,  data  were  taken  at  cooling  water  velocities 
given  by  20,  26,  35,  45,  54,  62,  70,  80,  and  20  percent  flow  rates.  For 
these  tubes,  the  80  percent  flow  rate  resulted  in  a  cooling  water 
velocity  of  about  4.4  m/s.  The  same  data  points  were  taken  for  the 
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medium-  and  large-diameter  tubes  using  steam  as  the  working  fluid. 
The  insert  in  the  small-diameter  tube  would  not  allow  flow  rates  above 
70  percent;  therefore,  for  this  tube  with  steam,  the  flow  rates  were 
adjusted  to  20.  26,  33,  40,  47,  54,  61,  66,  and  20  percent.  The  66 
percent  flow  rate  resulted  in  a  6.5  m/s  water  velocity  through  the 
tubes.  The  coolant  temperature  rise  decreases  with  increasing  water 
velocity,  thus  increasing  the  uncertainties  associated  with  this  tem¬ 
perature  measurement. 

Preliminary  tests  with  ethylene  glycol  indicated  that  sub-cooled 
boiling  was  occurring  inside  the  tube.  An  attempt  was  made  to  lower 
the  operating  pressure  of  the  system  to  decrease  the  vapor  tempera¬ 
ture.  However,  at  very  low  operating  pressures,  the  boiling  of  the 
ethylene  glycol  was  very  explosive  and  the  system  pressure  fluctua¬ 
tions  were  excessive  so  that  the  operation  was  considered  unsafe.  A 
combination  of  decreasing  operating  pressure  and  using  larger  flow 
rates  through  the  test  tube  was  used  to  prevent  the  internal  boiling. 
The  operating  pressure  was  lowered  to  62.2  mm  Hg.  which  resulted  in 
a  saturation  temperature  of  128°  C.  At  this  set  point,  the  system  still 
had  small  fluctuations  in  pressure  due  to  some  explosive  boiling.  The 
flow  rates  used  for  ethylene  glycol  were  30,  38,  45,  52,  60,  65,  70, 
80,  and  30  percent.  The  small-diameter  tubes  did  not  have  sufficient 
cooling  water  capacity  to  prevent  cooling  water  from  boiling  and 
therefore  were  not  tested.  Again,  the  coolant  temperature  rise 
decreases  with  increasing  water  velocity,  thus  increasing  the  uncer¬ 
tainties  associated  with  this  temperature  measurement.  Therefore, 
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owing  to  the  larger  velocities  involved,  the  data  for  ethylene  glycol  are 
expected  to  be  less  reliable  than  for  the  other  two  fluids. 

B.  DATA  REDUCTION  PROCEDURES 
1.  Background 

The  total  resistance  to  heat  transfer  from  the  vapor  to  the 
tube  cooling  water  consists  of  the  sum  of  the  vapor-side,  wall,  and 
inside  resistances,  as  shown  in  Figure  4.2. 


Tube 

Wall 


Wall 

Resistance 


Figure  4.1  Resistance  to  Heat  Transfer 


The  vapor-side  and  inside  resistances  were  convective  in  nature,  and 
therefore  the  resistances  could  be  expressed  by  the  reciprocal  of  the 
respective  heat-transfer  coefficient  and  surface  area  product. 


R,= 


1 

h|A,  * 


(4.1) 


and 


Rv  = 


1 


h0A0 


(4.2) 
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where 


a.  =  kD.(l  +  lit1i  +L2Tl2).  (4.3) 

and 

A0  =  xDrL  .  (4.4) 

Dr  =  root  diameter  (m) 

Ri  =  inside  resistance  (m2  •  K/W), 

Rv  =  vapor-side  resistance  (m2  •  K/W). 

hi  =  inside  heat- transfer  coefficient  (W/m2  •  K), 

h0  =  outside  heat-transfer  coefficient  (W/m2  •  K), 

At  =  effective  inside  area  (m2)1, 

A  =  outside  area  (m2), 

L  =  length  of  exposed  tube  (133  mm), 

Li  =  length  of  inlet  portion  of  tube  in  nylon  bushing  (m), 

L2  =  length  of  outlet  portion  of  tube  in  nylon  bushing  (m). 

Til  =  fin  efficiency  for  inlet  portion  of  tube, 
ri2  =  fin  efficiency  for  outlet  portion  of  tube. 

The  wall  resistance  was  easily  computed  from  the  thermal  conductiv¬ 
ity  of  copper  and  the  inside  and  outside  diameters. 

1The  axial  heat  conduction  into  the  inlet  and  outlet  insulated 
portions  of  the  tube  were  accounted  for  by  computing  the  associated 
fin  efficiencies.  For  this  purpose,  the  extended-fin  assumption  was 
used. 


lnT57 

Rw  -  Dr  ~v~: - 

w  r  2km 


(4.5) 


Rw  =  wall  resistance  (m2  •  K/W), 

km  =  thermal  conductivity  of  copper  (W/m  •  K). 

This  equation  assumes  that  the  heat  flux  was  radial  through  the  tube 
walls. 

Combining  Equations  4.1,  4.2,  and  4.5,  the  overall  thermal 
resistance  can  be  expressed  as: 

1 _ l _  i  Rw 

U0A0  "  h0A0+  hjA,+  A0  (4.6) 

In  order  to  compute  the  overall  heat-transfer  coefficient,  the 
heat- transfer  rate  must  be  found.  The  heat-transfer  rate  may  be  com¬ 
puted  by  measuring  the  inlet  and  outlet  test  tube  cooling  water  tem¬ 
peratures  and  computing  the  properties  of  the  cooling  water  at  its 
average  temperature.  The  test  tube  cooling  water  inlet  and  outlet  tem¬ 
peratures  were  measured  by  quartz  thermometers  and  a  ten-junction 
thermopile.  The  quartz  thermometer  readings  were  used  in  the  com¬ 
putation  of  the  log-mean-temperature  difference  (Equation  1.2)  and  in 
the  heat-transfer  rate  below: 


Q  =  mCpAT, 


(4.7) 
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where 


m=pcAcv,  (4.8) 

m  =  mass  flow  rate  (kg/s), 

Cp  =  specific  heat  at  constant  pressure  (J/kg  •  K), 

AT  =  cooling  water  temperature  rise  (K), 

Pc  =  test  tube  bulk  cooling  water  density  (kg/m3), 

Ac  =  cross-sectional  area  of  test  tube  (m2), 
v  =  test  tube  average  cooling  water  velocity  (m/s). 

Once  the  heat  transfer  rate  is  determined  by  Equation  4.7, 
the  value  can  be  substituted  into  Equation  1.2  to  determine  the  overall 
heat-transfer  coefficient.  The  determination  of  the  inside  and  outside 
coefficients  are  the  only  two  values  still  unknown  in  Equation  4.6. 

2.  Modified  Wilson  Plot  Technique 

In  general,  the  outside  resistance  should  be  minimized  if  the 
measurements  are  performed  only  for  the  inside.  Similarly,  the  inside 
resistance  should  be  minimized  if  the  primary  concern  is  the  outside 
heat-transfer  coefficient.  However,  the  modified  Wilson  plot  technique 
provides  both  inside  and  outside  coefficients  simultaneously.  With  this 
method,  in  order  to  maximize  accuracy,  the  inside  and  outside  resis¬ 
tances  should  be  made  as  equal  as  possible.  Therefore,  it  was  neces¬ 
sary  to  use  an  insert  to  boost  the  inside  coefficient  when  steam  was 
the  working  fluid. 

The  modified  Wilson  plot  technique  requires  that  the  form  of 
both  inside  and  outside  coefficients  be  known.  During  this  study. 


Sieder-Tate-type  and  Nusselt-type  equations  were  used  to  represent 
the  inside  and  outside,  respectively. 


(4.9) 


where 

kc  =  thermal  conductivity  of  cooling  water  (W/m  •  K), 

Ci  =  Sieder-Tate-type  coefficient. 

Re  =  Reynolds  number, 

Pr  =  Prandtl  number. 

Pc  =  dynamic  viscosity  of  cooling  water  at  bulk  temperature 
(N  •  s/m2), 

|iw  =  dynamic  viscosity  cooling  water  at  inner  wall  temperature 
(N  •  s/m2), 

and 


h0 


oi 


P,D0ATf 


tO  .25 


=  oF, 


(4.10) 


where 

a  =  dimensionless  coefficient, 

hfg  =  specific  enthalpy  of  vaporization  (J/kg), 

ATfg  =  average  temperature  difference  across  condensate  film  (K), 
Pi  =  dynamic  viscosity  of  condensate  (N  •  s/m2). 

By  substituting  the  Nusselt-  and  the  Sieder-Tate-type  equa¬ 
tions  into  Equation  4.6  and  rearranging,  the  following  is  derived: 
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[u  ■R-]F  = 


ACF  i 
C.ftA,  +  a 


(4.11) 


This  equation  is  in  linear  form  and  has  two  unknowns:  Ci  and  a.  By 


allowing: 


(4.12) 


a0f 

A,ft 


(4.13) 


C  — — 

i  “  m 


(4.14) 


(4.15) 


Equation  4. 1 1  reduces  to: 


Y  =  mX  +  b. 


(4.16) 


Even  though  it  is  not  readily  evident,  the  calculation  of  ft  and 
F  requires  that  Ci  and  a  be  known  because  of  fluid  property  variations 
with  temperature.  Therefore,  an  iterative  process  was  necessary  to 
compute  these  values.  To  begin  this  process,  a  value  of  2.5  was 
assumed  for  a  and  values  of  0.071  and  0.028  were  used  for  Ci  with  and 
without  an  insert,  respectively.  New  values  for  Ci  and  a  were  calcu¬ 
lated  by  performing  a  least-squares  fit  using  Equation  4.16.  This 
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process  was  repeated  until  consecutive  values  of  both  Ci  and  a  agreed 
within  ±0.1  percent. 

Once  the  Sieder-Tate-type  coefficient  was  computed,  the 
inside  heat-transfer  coefficient  was  determined  using  Equation  4.9. 
With  the  inside  and  overall  heat-transfer  coefficients  computed,  the 
outside  coefficient  was  easily  determined. 

_1 _ 1_ 

h0  ~U0 

Notice  that  any  error  in  computing  the  inside  coefficient  will  be  car¬ 
ried  over  into  the  value  of  the  outside  coefficient. 

Following  the  Nusselt  theory,  the  experimental  data  can  be 
expressed  and  be  fitted  using  a  least-squares  analysis  of  the  data  with 
the  following  form: 


h«A, 


+  R. 


f4  1  71 


q  =  aAT" 


(4.18) 


The  heat  flux  can  also  be  written  as  hAT,  resulting  in  the  following: 


h  =  aATn_1 


(4.19) 


According  to  the  Nusselt  theory,  n  =  0.75.  Therefore,  n  was  set  to  0.75 
in  this  investigation  so  that  the  enhancement  ratio  (based  on  constant 
vapor- side  temperature  drop)  can  be  expressed  as: 


€ 


AT 


hf 

h, 


<*fFf 

«>F,  ’ 


(4.20) 
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where  the  f  and  s  refer  to  finned  and  smooth  tubes,  respectively.  By 
keeping  the  average  temperature  drop  across  the  condensate  film  the 
same  for  the  finned  and  smooth  tube,  the  values  of  Ff  and  Fs  are  equal. 
Therefore.  Equation  4.20  becomes: 


€ 


AT 


(4.21) 


To  obtain  the  enhancement  ratio  at  constant  heat  flux,  the 
following  development  was  used: 


q  =  hf  AT  f  =  hs ATS 


(4.22) 


Therefore. 


e  q  - 


(4.23) 


From  Equation  4.18  with  n  =  0.75 


af  AT 


0.75 

f 


.  0.75 

asATs 


therefore 


(4.24) 


(4.25) 


and 


(4.26) 
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From  Equations  4.20  and  4.26,  it  is  clear  that 


-  r  l1  33 

€  q  _  l€  atJ  (4.27) 

Notice  that  with  the  above  development,  both  €q  and  cat  are  inde¬ 
pendent  of  q  and  AT. 

As  shown  by  Equation  4.10,  the  effect  of  tube  diameter 
appears  in  the  value  of  F,  and  therefore  a  is  independent  of  the  tube 
diameter  and  theoretically  should  not  vary.  Equation  4.11  shows  a  as 
the  inverse  of  the  y-intercept  with  values  of  X  and  Y  dependent  on  F. 
Experimentally  determined  values  of  a  for  smooth  tubes  may  differ 
slightly  from  one  tube  to  another.  However,  only  one  smooth  tube 
(having  a  diameter  of  19.05  mm,  Table  3.1)  was  available  during  this 
study  and  the  a  found  for  this  tube  (tube  #1)  was  used  in  computing 
enhancements. 

C.  MODIFICATIONS  TO  PROGRAM 

The  program  used  by  Zebrowski  [Ref.  12J  was  modified  for  this 
investigation.  The  modifications  incorporated  are: 

1.  The  tube  diameters  (both  inside  and  outside)  were  made 
variables. 

2.  Functions  that  compute  fluid  properties  were  extended  to 
include  ethylene  glycol. 

3.  Calibrations  were  included  to  account  for  the  viscous  dissipation 
from  the  insert  and  the  mixing  chamber  for  the  small-diameter 
tubes. 

The  ethylene  glycol  used  in  this  study  was  99  percent  pure  with  a 
range  of  194°  to  200°  C  for  the  normal  boiling  point.  Therefore,  it  was 
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not  possible  to  compute  the  non-condensing  gas  concentrations  accu¬ 
rately  from  the  experiential  measurements  described  previously. 
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V.  RESULTS  AND  DISCUSSION 


A.  INTRODUCTION 

Experimental  data  were  taken  as  described  in  Chapter  IV.  Some 
tubes  were  tested  at  least  twice  to  show  repeatability  of  the  apparatus 
used  during  this  investigation.  If  two  tests  did  not  result  in  outside 
heat-transfer  coefficients  within  ±5  percent,  a  third  set  of  data  was 
taken  on  the  same  tube.  Figures  5.1.  5.2,  5.3,  and  5.4  show  the 
repeatability  of  data  for  all  three  test  fluids.  As  shown  in  Figure  5.1, 
data  for  R-113  show  repeatability  within  ±3  percent  for  tubes  1  and 
18.  Similarly,  Figure  5.2  shows  repeatability  for  tube  1 1  with  steam  to 
be  within  ±7  percent.  This  repeatability  was  typical  of  all  the  medium- 
and  large-diameter  tubes  with  steam.  However,  poor  repeatability  with 
the  small-diameter  tubes  (tube  5)  with  steam  is  shown  in  Figure  5.3. 
The  variation  is  significantly  larger  at  ±12  percent.  As  seen  in  Figure 
5.4,  two  different  runs  of  ethylene  glycol  with  tube  19  yielded  con¬ 
densing  heat-transfer  coefficients  which  agree  within  ±5  percent. 

Owing  to  time  constraints,  a  method  of  recognizing  good  results 
without  necessarily  repeating  all  data  runs  needed  to  be  devised.  When 
taking  data,  the  flow  rate  in  the  test  tube  was  gradually  increased  from 
a  minimum  to  a  maximum  percent  and  a  last  set  of  data  was  taken 
again  at  the  minimum  flow  rate.  It  was  felt  that  a  comparison  of  the 
heat-transfer  coefficient  for  the  first  and  last  data  sets  was  a  good 
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Figure  5.1  Repeatability  of  R-113  Heat-Transfer  Coefficients 

for  Tubes  1  and  18 


Working  Fluid:  Ethylene  Glycol 

Pressure:  62  mm  Hg 


x 

\ 


CH  /°M 

Figure  5.4  Repeatability  of  Ethylene  Glycol  Heat-Transfer 

Coefficients  for  Tube  19 
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indication  of  the  quality  of  the  data.  If  the  comparison  showed  an 
increase  in  the  coefficient  from  start  to  finish,  dropwise  condensation 
could  have  been  responsible.  Similarly,  a  decrease  in  the  coefficient 
could  have  been  caused  by  an  in-leak  of  non- condensing  gases  (when 
operating  at  low  pressure),  and  by  the  generation  of  gases  within  the 
boiler  owing  to  the  decomposition  of  the  working  fluid.  As  stated  in 
Chapter  IV,  the  modified  Wilson  plot  technique  computed  the  Sieder- 
Tate-type  coefficient  by  using  a  linear  least- squares  fit  to  the  data  and 
taking  the  reciprocal  of  the  slope.  If  the  plot  of  the  data  did  not 
resemble  a  good  straight  line  (Figure  5.5),  the  data  run  was  discarded. 

B.  INSIDE  HEAT-TRANSFER  COEFFICIENT 

Table  5.1  lists  the  Sieder-Tate-type  coefficients  (Ci)  computed  for 
all  tubes  with  each  of  the  working  fluids.  As  discussed  earlier  in  Chap¬ 
ter  IV,  large-  and  medium-diameter  tube  families  shared  the  same 
insert,  while  a  different  insert  was  used  for  the  small  tubes.  Further, 
inserts  were  used  only  when  steam  was  the  working  fluid.  For  this 
reason,  Ci  values  were  computed  separately  for  the  small-diameter 
tubes.  As  can  be  seen,  the  modified  Wilson  plot  technique  gives 
slightly  different  Ci  values  for  different  tubes.  For  example.for  a  given 
family  of  tubes,  the  inside  diameter  was  identical  and  Ci  should  be  the 
same.  However,  as  shown  in  this  table,  Ci  for  steam  (low  pressure) 
varies  between  0.048  and  0.051.  This  variation  of  six  percent  was 
within  the  expected  uncertainty  from  the  modified  Wilson  plot  tech¬ 
nique.  Since  the  large  and  the  medium  tubes  have  identical  internal 
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005  0.01  0.015 


TABLE  5.1 


SUMMARY  OF  SIEDER-TATE-TYPE  COEFFICIENTS 


Tube 

s  (mm) 

R-113 

LP  Steam 

Working  Fluid 
Atm  Steam 

Ethylene 

Glycol 

Smooth  tube 

1 

- 

0.0278 

0.0679 

0.0663 

0.0318 

Small  tubes 

2 

0.25 

0.0279 

0.0494 

0.0509 

— 

3 

0.5 

0.0269 

0.0505 

0.0575 

— 

4 

1.0 

0.0281 

0.0507 

0.0513 

— 

5 

1.5 

0.0309 

0.0481 

0.0485 

— 

6 

2.0 

0.0395 

0.0505 

0.0482 

— 

7 

4.0 

0.0309 

0.0505 

0.0513 

— 

Medium  tubes 

8 

0.25 

0.0341 

0.0695 

0.0689 

0.0394 

9 

0.5 

0.0258 

0.0671 

0.0642 

0.0431 

10 

1.0 

0.0286 

- 

- 

0.0285 

11 

1.5 

0.0282 

0.0669 

0.0677 

0.0368 

12 

2.0 

0.0273 

0.0667 

0.0671 

0.0319 

13 

4.0 

0.0274 

0.091 

0.0671 

0.0299 

Large  tubes 

14 

0.25 

0.0359 

0.0626 

0.0668 

0.0425 

15 

0.5 

0.0279 

0.0632 

0.0624 

0.0350 

16 

1.0 

0.0317 

0.0636 

0.0684 

0.0336 

17 

1.5 

0.0303 

0.0631 

0.0687 

0.0327 

18 

2.0 

0.0301 

0.0666 

0.0679 

0.0318 

19 

4.0 

0.0302 

0.0668 

0.0678 

0.035 

QMC  tube 

20 

1.0 

.0293 

- 

0.0513 

- 

Working  Fluid 

R-113 
LP  Steam 
Atm  Steam 
Ethylene  Glycol 


0.030  ±0.002 
0.051  ±0.003 
0.066  ±0.002 
0.035  ±0.005 


(All  tubes) 

(Small  tubes) 

(Medium  and  large  tubes) 
(Medium  and  large  tubes) 
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geometries,  theoretically  they  should  have  the  same  Ci  value.  In  order 
to  obtain  the  outside  heat-transfer  coefficient  in  a  consistent  manner, 
the  average  Ci  values  were  computed  and  are  shown  at  the  bottom  of 
Table  5.1. 

For  R-113,  the  average  Sieder-Tate-type  coefficient  was  computed 
to  be  0.030  ±0.002.  As  seen  in  Table  5.1,  for  steam  at  atmospheric  and 
low-pressure  conditions,  nearly  the  same  Sieder-Tate-type  coefficient 
was  obtained.  Coefficients  for  the  medium  and  large  tubes  gave  an 
average  value  of  0.066  ±0.002.  The  small-diameter  tubes  (with  a  dif¬ 
ferent  insert)  gave  an  average  value  of  0.051  ±0.003  for  the  Sieder- 
Tate-type  coefficient. 

With  ethylene  glycol  (no  insert),  the  large  and  medium  tubes  gave 
an  average  Sieder-Tate-type  coefficient  of  0.035  ±0.005.  Unfortunately, 
it  was  impossible  to  prevent  subcooled  boiling  inside  the  small  tubes 
with  ethylene  glycol.  Notice  that,  as  discussed  by  Masuda  and  Rose 
[Ref.  3],  the  extent  of  subcooled  boiling  changes  with  the  water  veloc¬ 
ity.  For  this  situation,  it  was  not  possible  to  represent  the  inside  heat- 
transfer  coefficient  successfully  by  a  single  function  as  needed  for  the 
modified  Wilson  plot  technique.  Therefore,  the  testing  of  small  tubes 
was  discontinued  for  this  working  fluid. 

C.  OUTSIDE  HEAT-TRANSFER  COEFFICIENT 

1.  Repeatability  with  Previous  NPS  Data 

Figure  5.6  shows  the  repeatability  of  the  present  outside  heat- 
transfer  coefficient  for  R-113  as  compared  with  those  of  Zebrowski 
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Working  Fluid:  R-113  Present  Research  -Solid  line 

Pressure:  763  mm  Hg  Zebrowski  -Dash  line 


T - 1 - 1 - T - 1 - 1 - 1  8 
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Figure  5.6  Comparison  of  the  Present  R-113  Data  and  Zebrowski 

(Ref.  12] 
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[Ref.  12].  This  figure  shows  excellent  repeatability,  within  ±3  percent 
for  five  of  the  medium  tubes  (tubes  1.  9.  10,  11.  12,  and  13).  Similarly, 
Figures  5.7  and  5.8  show  repeatability  for  steam  under  low-pressure 
and  atmospheric  conditions,  respectively,  with  those  of  Lester  [Ref. 
13].  It  is  evident  that  the  maximum  disagreement  in  the  outside  coef¬ 
ficient  is  about  ±6  percent.  Figures  5.9  and  5.10  demonstrate  repeata¬ 
bility  with  results  from  Georgiadis  [Ref.  5]  with  steam  as  the  working 
fluid.  This  comparison  shows  a  disagreement  of  up  to  ten  percent  in 
the  outside  heat-transfer  coefficient. 

Notice  that  on  these  figures,  and  all  other  similar  figures, 
curves  through  the  data  are  drawn.  These  curves  represent  a  least- 
squares  fit  according  to  the  equation: 

q  =  aAT° 75  (5.1) 

2.  Outside  Heat-Transfer  Coefficient  for  R-113 

Figures  5.11  through  5.13  show  the  variation  of  the  outside 
heat-transfer  coefficient  for  R-113  versus  the  vapor-side  temperature 
drop  for  the  small,  medium  and  large  tubes.  On  each  figure,  the  curves 
are  identified  by  the  tube  number  (see  Table  3.1  for  the  corresponding 
fin  spacings).  Notice  that  the  relative  position  of  each  tube  remains  the 
same,  regardless  of  the  tube  diameter,  with  the  exception  of  the  0.25 
mm  fin  spacing  tubes  (tubes  2  and  8).  As  shown  by  the  uncertainty 
bands  in  Figure  5.11,  the  uncertainty  of  the  data  decreases  with  the 
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Figure  5.7  Comparison  of  the  Present  Low-Pressure 
Steam  Data  and  Lester  [Ref.  13) 
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Figure  5.10  Comparison  of  the  Present  Atmospheric 
Steam  Data  and  Georgiadis  [Ref.  5] 
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Figure  5. 


1  R-113  Heat-Transfer  Coefficients 
for  Small-Diameter  Tubes 
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Figure  5.13  R-113  Heat-Transfer  Coefficients 
for  Large-Diameter  Tubes 


increasing  vapor  temperature  drop.  Increasing  the  cooling  water  flow 
rates  increases  the  temperature  drop,  and  thereby  diminishes  the 
uncertainties.  This  is  true  for  all  fluids  tested. 

The  data  shown  in  these  three  figures  are  repeated  in  Table 
5.2  in  the  form  of  the  vapor-side  enhancement  ratios,  gat-  For  com¬ 
parison  purposes,  this  table  also  lists  the  area  enhancement  and  the 
condensate  retention  angle  for  each  tube.  The  area  enhancement  is 
the  ratio  of  the  total  surface  area  of  the  finned  tube  to  the  smooth  tube. 
It  is  evident  that  the  heat-transfer  enhancement  ratio  generally 
increases  with  decreasing  fin  spacing  (except  for  the  0.25  mm  spac¬ 
ing).  This  trend  is  easily  explained  by  the  increasing  area  enhance¬ 
ment  with  decreasing  fin  spacing.  However,  as  discussed  in  Chapter  II, 
a  decreased  fin  spacing  also  has  a  deleterious  effect  owing  to  the 
increased  condensate  retention  angle  (see  Table  5.2).  These  retention 
angles  were  computed  from  Equation  2.6.  Information  given  in  Table 
5.2  is  also  shown  graphically  in  Figure  5.14.  This  figure  shows  that,  for 
a  given  fin  spacing,  the  large-diameter  tube  gave  the  best  performance 
while  the  medium-diameter  tube  gave  the  poorest,  with  the  exception 
of  the  0.25  fin  spacing.  This  figure  also  shows  the  uncertainty  bands 
for  enhancements  predicted  (see  Appendix  B  for  details)  for  typical  R- 
113  data.  As  seen  in  Figure  5.14,  the  enhancements  for  each  family  lie 
within  or  very  close  to  the  uncertainty  band.  Although  the  uncertainty 
casts  doubt  on  the  relative  position  of  the  enhancements  displayed  on 
Figure  5.14,  it  is  important  to  recall  the  excellent  repeatability  with 
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TABLE  5.2 


SUMMARY  OF  R-113 


Fin  Area  Retention  Heat  Transfer 

Spacing  Enhancement  Angle1  Enhancement 


Tubes 

(mm) 

(degrees) 

(e  at) 

Small  tubes 

2 

0.25 

2.85 

180 

4.5 

3 

0.5 

2.54 

180 

5.2 

4 

1.0 

2.15 

156 

4.8 

5 

1.5 

1.93 

106 

4.5 

6 

2.0 

1.77 

87 

4.0 

7 

4.0 

1.46 

59 

2.8 

Medium  Tubes 

8 

0.25 

2.77 

79 

4.7 

9 

0.5 

2.47 

54 

5.2 

10 

1.0 

2.10 

37 

4.4 

1 1 

1.5 

1.88 

30 

4.4 

12 

2.0 

1.74 

26 

3.7 

13 

4.0 

1.44 

18 

2.7 

Large  Tubes 

14 

0.25 

2.73 

69 

5.3 

15 

0.5 

2.44 

47 

5.6 

16 

1.0 

2.08 

33 

5.0 

17 

1.5 

1.86 

27 

4.9 

18 

2.0 

1.72 

23 

4.5 

19 

4.0 

1.43 

16 

2.9 

9MC  Tube 

20 

1.0 

3.47 

43 

6.3 

1  Calculated  using  Equation  2.6. 
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FIN  SPACING  (mm) 


the  R-113  data.  Notice  that  among  the  tubes  tested,  the  tube  diameter 
has  no  effect  on  the  optimum  fin  spacing,  which  is  between  0.25  mm 
and  0.50  mm,  but  it  does  influence  the  magnitude  of  the  enhancement 
ratio. 

For  the  small- diameter  tubes  (Figure  5.11  and  Table  5.2), 
tube  2  (0.25  mm  spacing)  was  outperformed  by  tube  4  (1.0  mm  spac¬ 
ing).  As  the  diameter  is  increased  from  one  tube  family  to  another,  the 
enhancement  for  the  0.25  mm  tube  overtakes  the  enhancement  of  the 
1.0  mm  tube.  In  the  medium-diameter  tubes,  the  0.25  mm  tube 
(tube8)  performs  slightly  better  than  the  1.0  mm  (tube  10).  This  trend 
continues  into  the  large-diameter  family,  where  tube  14  (0.25  mm 
spacing)  outperformed  tube  16  (1.0  mm  spacing)  by  a  larger  magni¬ 
tude.  As  seen  in  Table  5.2,  the  trend  is  consistent  with  the  decreasing 
condensate  retention  angle.  Comparing  the  retention  angles  for  the 
small  to  large  tubes  reveal  that  the  condensate  retention  angle 
decreases  by  30  degrees  for  the  0.25  mm  and  only  20  degrees  for  the 
1.0  mm  tube.  Therefore,  as  the  diameter  increases  from  small  to  large, 
the  unflooded  area  for  tubes  with  a  0.25  mm  fin  spacing  increases 
much  more  than  that  for  the  tubes  with  a  1.0  mm  fin  spacing. 

Tests  with  the  QMC  tube  resulted  in  an  enhancement  ratio 
(for  constant  vapor-side  temperature  drop)  of  6.3.  As  summarized 
graphically  by  Masuda  and  Rose  [Ref.  3),  the  1.0  mm  QMC  tube  as 
tested  by  Yau,  et  al.  [Ref.  4],  with  a  0.7  m/s  vapor  velocity,  had  approx¬ 
imately  an  enhancement  of  6.4.  Thus,  the  agreement  between  these 
two  independent  investigations  is  excellent. 
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3.  Outside  Heat-Transfer  Coefficient  for  Steam 


Figures  5.15  through  5.20  show  the  variation  of  the  outside 
heat-transfer  coefficient  for  steam  for  the  small-,  medium-,  and  large- 
diameter  tubes.  These  figures  were  constructed  in  the  same  manner 
as  those  for  R-113.  Notice,  as  with  R-113,  the  uncertainty  is  much 
greater  for  the  low  flow  rates  (Figure  5.18).  As  the  difference  between 
the  vapor  temperature  and  outside  wall  temperature  decreases,  the 
uncertainty  increases.  The  large-diameter  tubes  have  the  greatest 
uncertainties  (see  Appendix  B)  due  to  the  very  low  temperature 
difference.  Owing  to  this  low  temperature  difference  and  associated 
large  uncertainty,  the  reprocessed  data  for  the  large-diameter  tubes 
with  an  average  Cj,  listed  in  Table  5.1,  was  not  well -correlated  by 
Equation  5.1,  which  uses  an  exponent  of  0.75  on  the  vapor  to  wall 
temperature  difference.1  A  summary  of  the  vapor-side  enhancements, 
retention  angles,  and  the  area  enhancements  are  located  in  Tables  5.3 
and  5.4. 

Enhancement  versus  fin  spacing  for  steam  at  low-pressure 
conditions  is  shown  in  Figure  5.21.  Using  Equation  2.6,  the  fin  spacing 
at  which  each  tube  diameter  becomes  fully  flooded  was  computed  and 
indicated  on  the  figure  by  a  vertical  line.  Notice  that  the 


1  A  least  squares  fit  of  the  data  using  Equation  5.1  but  allowing  a 
different  exponent  yielded  exponents  of  0.45  and  0.61  for  the  large 
diameter  tubes  at  low-pressure  and  atmosphere  conditions, 
respectively. 
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Figure  5.15  Steam  Heat-Transfer  Coefficients  for 
Small-Diameter  Tubes  at  Low-Pressure  Conditions 
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Figure  5.18  Steam  Heat-Transfer  Coefficients  for 
Sm'U-Diaxneter  Tubes  at  Atmospheric  Conditions 
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TABLE  5.3 


SUMMARY  OF  STEAM  ENHANCEMENTS 
AT  LOW-PRESSURE  CONDITIONS 


Tubes 

Pin 

Spacing 

(mm) 

Area 

Enhancement 

Retention 

Angle1 

(degrees) 

Small  tubes 

2 

0.25 

2.85 

180 

3 

0.5 

2.54 

180 

4 

1.0 

2.15 

156 

5 

1.5 

1.93 

106 

6 

2.0 

1.77 

87 

7 

4.0 

1.46 

59 

Medium  Tubes 

8 

0.25 

2.77 

180 

9 

0.5 

2.47 

180 

10 

1.0 

2.10 

109 

1 1 

1.5 

1.88 

84 

12 

2.0 

1.74 

71 

13 

4.0 

1.44 

48 

Large  Tubes 

14 

0.25 

2.73 

180 

15 

0.5 

2.44 

180 

16 

1.0 

2.08 

92 

17 

1.5 

1.86 

72 

18 

2.0 

1.72 

61 

19 

4.0 

1.43 

42 

1  Calculated  using  Equation  2.6. 


Heat  Transfer 
Enhancement 

(sat) 


2.1 

2.0 

2.3 

2.0 

2.2 

2.1 


1.9 

1.8 

2.6 

2.5 

2.4 


2.2 

1.8 

2.5 

3.0 

2.5 

2.5 
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TABLE  5.4 


SUMMARY  OF  STEAM  ENHANCEMENTS 
AT  ATMOSPHERIC  CONDITIONS 


Tubes 

Fin 

Spacing 

(mm) 

Area 

Enhancement 

Retention 

Angle1 

(degrees) 

Heat  Transfer 
Enhancement 

(6  at) 

Small  tubes 

2 

0.25 

2.85 

180 

3.1 

3 

0.5 

2.54 

180 

2.8 

4 

1.0 

2.15 

134 

2.4 

5 

1.5 

1.93 

97 

2.3 

6 

2.0 

1.77 

81 

2.5 

7 

4.0 

1.46 

55 

2.2 

Medium  Tubes 

8 

0.25 

2.77 

180 

2.4 

9 

0.5 

2.47 

180 

2.2 

10 

1.0 

2.10 

100 

_ 

11 

1.5 

1.88 

78 

2.9 

12 

2.0 

1.74 

66 

2.7 

13 

4.0 

1.44 

45 

2.3 

Large  Tubes 

14 

0.25 

2.73 

180 

2.7 

15 

0.5 

2.44 

147 

2.4 

16 

1.0 

2.08 

85 

2.8 

17 

1.5 

1.86 

67 

3.1 

18 

2.0 

1.72 

57 

3.0 

19 

4.0 

1.43 

40 

2.9 

QMC  Tubes 

20 

1.0 

3.47 

124 

3.6 

1  Calculated  using  Equation  2.6. 
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Figure  5.21  Enhancements  for  Steam  at  Low-Pressure  Conditions 
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small-diameter  tube  becomes  fully  flooded  at  a  spacing  just  under 
1.0  mm.  Medium-  and  large-diameter  tubes  become  fully  flooded  at  a 
spacing  of  approximately  0.5  mm.  In  the  fully  flooded  region, 
enhancement  was  expected  to  increase  as  the  fin  spacing  decreased. 
Once  a  tube  is  fully  flooded,  any  further  decrease  in  fin  spacing  should 
replace  the  flooding  condensate  with  tube  material  and  result  in  a 
greater  performance.  This  trend  is  seen  in  the  medium-  and  large- 
diameter  tubes.  Figure  5.21  indicates  that,  for  small-diameter  tubes, 
the  enhancement  in  the  fully  flooded  region  appears  to  be  more  com¬ 
plicated  than  simply  replacing  condensate  with  tube  material. 

For  partially  flooded  tubes,  the  change  in  diameter  from 
medium  to  large  did  not  affect  the  optimum  fin  spacing  (1.5  mm). 
However,  the  small -diameter  tube  exhibited  an  optimum  fin  spacing 
close  to  2.0  mm.  Furthermore,  for  partially  flooded  tubes.  Figure  5.21 
shows  that,  for  a  given  fin  spacing,  the  large-diameter  tube  gave  the 
best  performance  while  the  small-diameter  tube  gave  the  poorest 
performance. 

Figure  5.22  displays  the  enhancement  for  steam  at  atmo¬ 
spheric  conditions.  Again,  the  vertical  lines  indicate  the  fin  spacing  at 
which  the  tube  becomes  fully  flooded.  Large  increasing  enhancements 
are  shown  to  the  left  of  the  fully  flooded  lines  (decreasing  fin  spacing). 
As  with  the  low-pressure  conditions,  the  optimum  fin  spacing  is 
1.5  nan  for  the  medium-  and  large-diameter  tubes,  while  a  2.0  mm 
optimum  fin  spacing  is  shown  for  the  small-diameter  tube.  The 
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Figure  5.22  Enhancement*  for  8team  at  Atmospheric  Condition* 
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optimum  fin  spacing  for  the  medium -diameter  tube,  in  both  low- 
pressure  and  atmospheric  conditions,  agrees  with  Georgiadis.  The 
order  of  enhancement  from  least  to  greatest  is  small,  medium,  and 
large. 

The  QMC  tube  is  shown  in  Figure  5.22  to  have  an  enhance¬ 
ment  of  3.5.  This  is  far  larger  than  the  measured  enhancement  of  only 
2.2  obtained  at  Queen  Mary  College  [Ref.  3].  Marto  [Ref.  30]  pointed 
out  in  his  review  paper  that  there  was  an  unexpected  inconsistency 
between  the  NPS  obtained  data  and  the  Queen  Mary  College  obtained 
data  for  R-113  and  steam.  He  found  that  for  R-113,  Queen  Mary  Col¬ 
lege  data  [Ref.  3]  were  higher  than  the  NPS  data  (Zebrowski  [Ref.  12]). 
As  pointed  out  earlier  in  this  thesis,  for  R-113,  the  present  results 
with  the  QMC  tube  (this  investigation)  are  in  excellent  agreement  with 
the  data  obtained  at  Queen  Mary  College.  Marto  (Ref.  301  also  pointed 
out  that  for  the  steam  at  atmospheric  pressure,  the  data  taken  at 
Queen  Mary  College  were  less  than  the  data  taken  a  NPS.  However,  for 
steam,  the  result  of  this  thesis  for  the  enhancement  of  the  QMC  tube 
was  much  larger  than  the  data  for  the  NPS  tubes.  This  indicates  a  con¬ 
sistent  trend  for  R-113  and  steam  when  comparing  the  performance 
of  tubes  4  and  20.  The  reasons  for  the  inconsistent  trend  of  the  Queen 
Mary  College  data  are  not  known  at  present.  Further  investigation  and 
communication  is  needed  to  solve  this  discrepancy. 
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4.  Outside  Heat-Transfer  Coefficient  for  Ethylene  Glycol 

Figures  5.23  through  5.24  display  the  variation  of  the  outside 
heat-transfer  coefficient  with  the  vapor-side  temperature  drop  for  the 
medium-  and  large-diameter  tubes.  These  figures  were  constructed  in 
the  same  manner  as  the  R-113  figures.  The  relative  position  of  each 
tube  remained  unchanged  for  each  diameter.  The  uncertainties  for  the 
ethylene  glycol  are  associated  with  the  standard  deviation  of  experi¬ 
mental  values  of  Cj.  Even  though  the  uncertainty  of  the  ethylene  glycol 
data  is  as  high  as  the  large-diameter  steam  data.  Figures  5.23  and  5.24 
show  good  agreement  with  the  least  squares  fit  of  Equation  5.1.  As 
shown  in  Figures  5.23  and  5.24.  the  medium  tube  slightly  outper¬ 
formed  the  large  tube.  This  is  contradictory  to  the  trends  of  R-113 
and  steam.  As  shown  in  Figure  5.25.  the  difference  in  enhancement 
between  the  medium-  and  large-diameter  tubes  is  within  the  uncer¬ 
tainty  band,  and.  therefore,  the  large-diameter  tubes  could  possibly 
have  the  better  enhancement.  In  either  case,  the  performance  of  the 
two  tube  families  is  approximately  the  same.  A  summary  of  the  vapor- 
side  enhancements,  retention  angles,  and  area  enhancements  is 
located  in  Table  5.5.  The  change  in  diameter  from  medium  to  large 
did  not  affect  the  optimum  fin  spacing  of  1.0  mm  for  the  tubes  tested. 
This  optimum  fin  spacing  agrees  with  Masuda  and  Rose  [Ref.  31. 

D.  SUMMARY  OF  RESULTS 

The  literature  discussed  in  Chapter  II  indicated  that  an  increase 
in  root  diameter  will  decrease  the  condensate  retention  angle,  and. 
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>  /K 


Medium  Tubes 
Larqe  Tubes 


ENHANCEMENTS  (AT) 


Figure  5.25  Enhancements  for  Ethylene  Glycol 
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FIN  SPACING  (mm) 


TABLE  5.5 

SUBSMART  OP  ETHYLENE  GLYCOL  ENHANCEMENTS 


Fin 

Area 

Retention 

Heat  Tranafer 

Spacing 

Enhancement 

Angie1 

Enhancement 

Tubes 

(mm) 

(degrees) 

(gat) 

Medium  Tubes 

8 

0.25 

2.77 

180 

1.3 

9 

0.5 

2.47 

102 

2.3 

10 

1.0 

2.10 

66 

2.6 

1 1 

1.5 

1.88 

53 

2.3 

12 

2.0 

1.74 

46 

2.1 

13 

4.0 

1.44 

32 

1.6 

Large  Tubes 

14 

0.25 

2.73 

151 

1.5 

15 

0.5 

2.44 

86 

2.2 

16 

1.0 

2.08 

58 

2.5 

17 

1.5 

1.86 

47 

2.3 

18 

2.0 

1.72 

40 

2.0 

19 

4.0 

1.43 

28 

1.7 

1  Calculated  using  Equation  2.6. 
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therefore,  a  clear  trend  of  increasing  enhancement  from  small-  to 
large-diameter  tubes  was  expected.  R-113  had  similar  enhancements 
for  the  small-  and  medium-diameter  tubes,  while  the  ethylene  glycol 
had  similar  enhancements  for  the  medium  and  large.  The  possibility  of 
competing  effects  exists  with  a  change  in  root  diameter.  For  the  large- 
diameter  tubes,  condensate  must  flow  along  a  longer  path  length  from 
the  top  to  the  bottom  of  the  tube.  This  longer  path  length  yields  a 
larger  average  film  thickness  in  the  unflooded  portion  of  the  tube 
when  compared  with  a  small  tube.  This  degrades  the  performance  in 
the  unflooded  portion  of  the  large-diameter  tube  when  compared  to  a 
small-diameter  tube.  The  significance  of  this  effect  depends  on  fluid 
properties  such  as  thermal  conductivity,  viscosity,  and  surface  tension 
and  warrants  further  study. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

1.  Within  the  range  of  diameters  tested,  the  effect  of  root  diameter 
on  the  vapor-  side  coefficient  was  small.  Results  indicate  that  two 
or  more  competing  mechanisms  may  exist  while  increasing  the 
root  diameter. 

2  With  the  exception  of  the  small -diameter  tubes  with  steam,  the 
optimum  ftn  spacing  was  near  0.5.  1.0.  and  1.5  mm  for  R-l  13. 
ethylene  glycol,  and  steam  respectively. 

2  The  optimum  ftn  spacing  for  the  small -diameter  tube  with  steam 
was  approximately  2.0  mm. 

4.  With  the  exception  of  the  small -diameter  tube  with  steam,  a 
change  in  diameter  did  not  effect  the  optimum  fin  spacing  for 
each  test  fluid. 

S  The  data  presented  with  this  Investigation  supports  the  findings 
of  oast  research  at  NFS.  Variations  from  previous  data  for  R-l  13 
and  steam  were  within  ±3  percent  and  ±10  percent  respectively. 

E  The  comparison  of  the  enhancements  for  the  QMC  tube  was  out¬ 
standing  for  R-l  13.  while  the  steam  comparison  exhibited  a  large 
discrepancy. 


B.  RECOMMENDATIONS 

1.  Retest  the  families  of  tubes  to  verify  trends  observed  for  ethylene 
glycol. 

2  Manufacture  a  new  tube  with  the  same  dimensions  as  tube  10 
using  the  same  company  that  manufactured  the  original  tubes. 

3.  Manufacture  smooth  tubes  of  small  and  large  diameters.  Compare 
the  values  obtained  of  a  for  these  small  and  large  smooth  tubes 
with  the  value  of  a  for  the  medium  tube. 

4.  Manufacture  additional  tubes  to  include  more  fin  spacings.  espe¬ 
cially  in  the  vicinity  of  the  fully  flooded  condition. 
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5.  Operate  the  apparatus  with  methanol  or  ethanol  as  working  flu¬ 
ids.  In  addition  to  having  low  boiling  temperatures  and  good  wet¬ 
ting  characteristics,  these  fluids  have  approximately  the  same 
surface-tension  to  density  ratio  as  ethylene  glycol. 

6.  Increase  the  cooling  water  flow  rates  through  the  test  tube  to 
minimize  uncertainties  in  vapor-side  heat-transfer  coefficients. 

7.  Modify  the  apparatus  with  a  controller  for  the  auxiliary  condenser 
cooling  water.  Once  the  equilibrium  set  temperature  is  obtained, 
the  controller  can  be  placed  in  operation  to  maintain  the  set 
temperature. 

8  Provide  a  cooling  water  sump  and  pump  for  the  auxiliary  con¬ 
denser.  The  present  tap  water  system  is  very  susceptible  to  pres¬ 
sure  fluctuations  in  the  building. 


APPENDIX  A 

RAW  DATA 


Table  A.  1  contains  names  of  raw  data  files  with  corresponding 
tube  number.  Actual  raw  data  files  follow  table. 
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21.90 

48.32 

C 

1.9? 

20.79 

2!  .51 

48. 3E 

6 

I  .  9? 

20.79 

21.51 

49.34 

7 

:.si 

20.66 

21  .30 

48.46 

6 

2.5' 

20. G9 

21.30 

48.45 

5 

3.0C 

20.62 

21.16 

49.47 

*i  c 

3 . 0i3 

20.62 

„  '  .  1  ' 

40.46 

i  ' 

3.43 

20.55 

21  .0? 

48.44 

i  2 

3.43 

20.59 

2l.CS 

46 . 44 

•  *? 

7  i  " 

20.ee 

2 1  . 00 

46.37 

1  4 

3 . 5“ 

20.55 

2'  .00 

4  8.37 

i  = 

4 . 4  '• 

“*  g  C  1 

20.52 

4E  .44 

1  £ 

4 . 4  > 

20.5' 

20 . 52 

48.44 

•  “* 

*  .  1  E 

2  '  ,  2 

21 . 2~ 

46 .44 

1  £ 

•  '  £ 

2  ■  .2" 

n  **  "i  — 

48.44 

File 

« •*  ♦  : 

L5F-10 

» 

c'es 

a  r*  •  C  Z  r  3 

i  *.  1  or  : 

-t wispier  i 

c 

v  a:  : 

r*  l_  *  ’  -  -  j  • 

.• : 

.  2r  *  5 

« 

2  *  ♦ » 

T  l  r 

TPut 

T  t 

t 

'  a 

-  ' 

;  - 

'  •-  • 

1  .  ’  E 

2  '•  .  ’  E 

ee 

46.41 

’  .  *  t 

2  ‘  .  •  = 

2’  .96 

46 .44 

z 

1.4  = 

20.52 

2'  .92 

4S  .  44 

4 

i  .4  = 

20.51 

2’  .61 

4S .  46 

I  .  5*" 

ze.7e 

21.27 

46.36 

6 

i  .9? 

20.70 

21.27 

48.34 

7 

2.51 

20.59 

21.06 

48.40 

c 

2.51 

20.56 

21.06 

48.35 

s 

3.00 

20.53 

20.95 

48.35 

tc 

3.00 

20.53 

20.95 

48.36 

1 1 

3.43 

20.50 

20.87 

48.41 

i : 

3 .43 

20.50 

20.68 

48.42 

!  3 

3. 87 

20.45 

20.93 

48.48 

1  4 

3.6r 

20.45 

20.93 

46.51 

1  5 

4.4’ 

20. 4£ 

20.76 

48.45 

’E 

4.4’ 

20. 4£ 

20.  ”E 

46.48 

•  - 

;  .  ’  £ 

,  -*  p 

k  »  .  *.  «• 

22.05 

48.40 

• 

1 8 

“  •  m  0 
«* 

22.06 

46.40 
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File  Nane  :  Ml AEGV4S 

Pressure  Condition:  Vacuur 
Vapor  Velocity:  10.0  (n/s) 


Data 

Vw 

Tin 

Tout 

Ts 

* 

<  1*1/  s  ) 

1  C  > 

i  C  ) 

(C  ) 

1 

1  .70 

22.12 

24.63 

127.86 

~y 

1  .70 

22.12 

24.65 

127.68 

3 

z.  13 

22.00 

24.06 

127.92 

4 

2.13 

22.00 

24.07 

127.91 

5 

Z.E1 

21.92 

23.71 

127.93 

5 

2.51 

21.92 

23.71 

127.90 

“7 

2.89 

21  .87 

23.44 

127.87 

8 

2.89 

21.87 

23.45 

127.88 

Q 

o .  -  2 

21.82 

23.21 

127.97 

'0 

0.0^ 

21.82 

23.21 

127.96 

1  1 

3 .59 

^  1  “IQ 

L.  1  .  (O 

23.03 

127.94 

i ; 

3.55 

2  t  #  75 

23.09 

127.96 

i  ™ 

3.86 

21 

22.98 

127.93 

1  4 

3.6c 

01.?? 

22.98 

127.91 

1  5 

4  .  40 

21.73 

22.79 

127.51 

16 

4 . 40 

21.73 

ya 

127.90 

1 ~ 

1 .70 

“i*:  *  "t 

•  <  • 

24.  ES 

127.96 

^  c 

1  .70 

->  j  *7 

24.58 

127.54 

File 

Name : 

M4E5V39 

Pres 

sure  Ccnc 

,  4 

1  C  r  : 

Vacuur. 

Vapc 

-  veiocit 

>  : 

10.0  (  1*1  / 

5  > 

Cat  s 

Vw 

i  :  n 

Tout 

Ts 

£ 

nr  C 

'  C  > 

•  2  '■ 

f  c 

1 

Z  .  I 

3 

i  .36 

24.52 

1 2". 91 

Z 

-  r 

•  .36 

24.51 

127.95 

7 

Z  .  ST 

'  .  3! 

24.11 

127.93 

4 

Z  .  5” 

’  .31 

24.12 

127.9! 

c 

Z .  £5 

1  .25 

23.60 

527.93 

6 

2.85 

1.25 

23.80 

127.83 

7 

3.21 

*» 

1.20 

23.52 

127.92 

8 

3.21 

** 

4. 

i.20 

23.50 

127.97 

9 

3.59 

n 

1.16 

23.26 

127.92 

10 

3.59 

n 

i.15 

23.25 

127.99 

1  1 

3.57 

z 

1.12 

23. C4 

127.93 

1  2 

3.97 

z 

1.12 

23.05 

128.00 

i  7 

4.30 

z 

1  .09 

22.88 

127.91 

14 

4 . 30 

4_ 

1  .09 

22.86 

127.87 

15 

4.5“ 

n 

1  .05 

22.71 

127.91 

’  £ 

4.6“ 

n 

i  .05 

22. "0 

127.93 

;  ? 

2 . 24 

i  .34 

24.46 

127.93 

1  5 

2.24 

i  .  34 

24.44 

127.68 
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F 1  ie 

Marne : 

M5ESV34 

Press 

ure  Condi 

t  i  o  r, : 

Vacuum 

Vapor 

Veloc 1 t  y 

; 

10.0  ( m 1 s 

; 

Data 

Vw 

Tin 

Tout 

Ts 

t  '' 

(  M  /  5  ) 

(C  ) 

(C  ) 

(C  ) 

1 

2.79 

20.45 

23.23 

127.79 

-> 

2.79 

20.45 

23.22 

127.79 

3 

3.11 

20.40 

22.96 

127.78 

4 

3.11 

20.40 

22.97 

127.84 

5 

3.44 

20.36 

22.75 

127.81 

6 

3.44 

20.36 

22.74 

127.85 

7 

3 . 75 

20.33 

22.55 

127.60 

8 

3.76 

20.33 

22.55 

127.79 

9 

4.19 

20.29 

22.34 

127.80 

a? 

4.19 

20.29 

22.34 

127.78 

!  ! 

4.53 

20.26 

22.15 

127.79 

1  z 

4.53 

20.26 

22.15 

127.82 

1  w’ 

5.05 

20.23 

21 .99 

127.84 

'  4 

5.05 

20.23 

21.99 

127.82 

’5 

5 . 3S 

20.21 

21  .8" 

127.80 

IE 

5.38 

20.21 

21  .59 

127.81 

17 

2.79 

20 . 45 

23.24 

127.75 

1  8 

2.79 

20.45 

23.24 

1 27 . 79 

File  Mens:  M6EGV40 

Condition:  Vacuum 
vapor  Velocity:  1  0 .  0  ( m  s 


Data 

Vw 

T 

l  r. 

T  out 

Ts 

t 

:  f  5  ■’ 

< 

C  > 

(  C  1 

(  C  > 

1 

J.  • 4 

2  ! 

.42 

24.55 

127.90 

J. 

2.24 

21 

.43 

24.56 

127. 92 

2.57 

•’  i 

.39 

24.18 

127.92 

4 

2.57 

2  1 

.  40 

24.16 

127.97 

5 

2.89 

2! 

.36 

23.89 

127.87 

6 

2.89 

21 

.36 

23.89 

127.86 

n 

3.21 

21 

.34 

23.66 

127.90 

8 

3.21 

21 

.34 

23.68 

127.91 

9 

3.59 

21 

.31 

23.43 

127.89 

10 

3.59 

21 

.31 

23 . 4o> 

127.91 

1  1 

3.97 

21 

.29 

23.24 

127.85 

1  2 

3.97 

21 

.29 

23.24 

127.84 

13 

4.29 

21 

.28 

23.10 

127.91 

1  4 

4.29 

21 

.28 

23.10 

127.95 

15 

4.67 

21 

.26 

22.96 

127.91 

16 

4.57 

2  i 

.26 

22 . 9E 

127.86 

17 

2.24 

2  i 

.57 

24.65 

127.92 

1  5 

2 . 24 

2  i 

■ 3  ' 

24.70 

126 

127.89 

File 

Name : 

M7EGV41 

Pressure  Condit 

ion: 

Vacuum 

Vapor 

Veloci 

10.0  ( m : 5 

) 

Data 

Vw 

Tin 

Tout 

Ts 

« 

(  m ,  s  ) 

(C  > 

<c; 

<  C  ) 

1 

2.24 

n 

1  .51 

24.35 

128.00 

n 

2.24 

n 

1 .51 

24.36 

127.97 

3 

2.57 

4» 

1  .45 

24.00 

127.97 

4 

2.57 

n 

L 

1  .45 

24.01 

128.03 

5 

2.85 

n 

1 .40 

23.72 

127.96 

6 

2.89 

n 

1  .40 

23.73 

127.99 

7 

3.21 

n 

1  .35 

23.49 

127.96 

6 

3.21 

*■> 

1  .35 

23.50 

127.99 

9 

3.55 

i. 

1  .31 

23.28 

128.01 

10 

3.59 

1 .31 

4.  Z)  .  4-  ( 

127.99 

’  t 

3.97 

2 

1  .28 

23.06 

127. 98 

1  2 

3.9n 

n 

1  .28 

23.08 

t  27 . 95 

’3 

4.29 

2 

1  .25 

22.92 

127.93 

1  4 

4.29 

2 

1  .25 

22 . 92 

127.97 

i  c 

4. El1 

n 

1  2  2 

22.75 

128.00 

16 

4 .  E7 

n 

1  .22 

22.79 

127.95 

!  7 

2.24 

n 

1  .52 

24.39 

127.97 

1  c 

2.24 

2 

1  .52 

24.39 

127.97 

File 

Name : 

MSEGV4Z 

Pres 

sure  Condi 

ti 

or, : 

Vacuum, 

Vapo 

r  Velocity 

| 

10.0  U»'5 

Data 

V* 

T 

1  r 

T  out 

Ts 

s 

1  m ,  s  ) 

i 

c 

l  C  ) 

(  C  ; 

1 

2.24 

2 1 

.  46 

23.73 

127.97 

*■* 

“  n  a 

2  1 

.50 

L  4  .  S  *. 

127.98 

■7 

2 . 5" 

21 

.47 

23.53 

127.92 

4 

2.5"’ 

2  1 

.  47 

23.53 

127.91 

5 

2.85 

21 

.43 

23.30 

127.92 

6 

2.35 

21 

.  43 

23.31 

127.94 

7 

3.21 

21 

.40 

23.12 

127.99 

8 

3.21 

21 

.40 

23.  1  1 

127.96 

9 

3.59 

21 

.36 

22.93 

128.02 

10 

3.59 

21 

.36 

22.92 

128.00 

1  1 

3.97 

21 

.  33 

22.77 

127.94 

12 

3.97 

21 

.33 

22.75 

127.90 

13 

4.29 

2-1 

.31 

22.65 

127.91 

1  4 

4.29 

21 

.31 

22.65 

127.93 

15 

4.6? 

2  1 

.25 

22.54 

127.93 

16 

4.67 

.  29 

22.53 

127.93 

1  *“ 

2 . 24 

w  i 

_  C  q 

2  3.95 

>27.58 

{  C 

2 . 24 

I  1 

.60 

23.95 
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1 Z? . 5C 

File  Name : 

M 1 EGV37 

F'ressu 

re  Condi 

t  ion: 

Vacuum 

Vapor 

Veloc i t y 

I 

10.0  ( m / 

Data 

Vw 

Tin 

Tout 

t 

(  M  /  S  i 

(C) 

tC  ) 

1 

1  .  16 

7 9 

43 

25.23 

2 

1.18 

9  2 . 

45 

25.27 

3 

1  .48 

9  9 

28 

24.56 

4 

1.43 

n 

im  t- 

28 

24.57 

5 

1  .97 

9  n 

*2  L- 

12 

23.91 

G 

1  .97 

n 

«2  L. 

12 

23.91 

7 

2.51 

9  9 

.01 

23.45 

8 

2.51 

9  *"1 

.02 

23.47 

9 

3.00 

21 

.95 

23.  IS 

10 

3.00 

2 1 

.95 

23.19 

1  1 

3.43 

21 

.91 

23.00 

1  2 

3.43 

2  1 

.91 

23.00 

1  3 

3.85 

21 

.87 

22.85 

1  4 

3 .  £8 

21 

.8? 

22.85 

15 

4.40 

9  1 

.  84 

22.69 

1  5 

4.40 

2 1 

.84 

22.70 

i  7 

1.16 

2  *2 

.57 

25.41 

1  8 

1.16 

L. 

.57 

25.40 

Ts 

(C  > 

127.78 
127.63 
127.85 

127.83 

127.79 
127.78 

127.82 

127.83 

127.84 

127.82 

127.82 

127.75 
127.63 

127.82 
127.78 
127.81 
127. 7G 

127.76 
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ai  m  i_.j  t  )  'Smon-sioimt'  t.i  t  i  «  a>  w  i  m  ~i  cn  in  t>  w  r  >  —  1S1  ud 


File 

Name  : 

L1EG53 

Press 

u-e  Coro 

1 1 1 or  : 

Vacuum 

Vapor 

Veioc : i 

y  • 

10.0  (1*1/5 

Dat  a 

VlAi 

Tin 

Tout 

Ts 

* 

(1*1/5  .’ 

(C  ) 

(C  ) 

(C  ) 

1 

1  .70 

nn  -in 
ti  .  (  L. 

25.76 

127.99 

2 

1  .70 

22.71 

25.77 

127.39 

3 

2.13 

22.57 

25.10 

127.93 

4 

2.13 

22.58 

25.09 

127.91 

5 

2.51 

22.48 

24.67 

127.98 

g 

2.51 

22.48 

24.68 

128.05 

7 

2.89 

22.41 

24.35 

127.98 

8 

2 .  £9 

22.4! 

24.36 

128.05 

9 

3 . 32 

22.35 

24.07 

127.99 

1  £■ 

3.32 

22.35 

24.08 

128.05 

}  \ 

3 . 59 

-i  -> 

23.93 

127.99 

12 

3 .53 

->  “  "7  n 

*7  o*  "7 
<u  _>  •  w 

127.96 

;  ^ 

5.  S3 

22.23 

23.79 

128.00 

1  4 

3 . 6E 

22.29 

23.79 

127.97 

1  c 

4.40 

22 . 25 

23.57 

128.01 

1  b 

4 . 40 

22.25 

23.57 

127.91 

*  *7 

1  .  70 

22.70 

25.79 

128.03 

1  8 

1  .70 

22.70 

25.79 

127.97 

1 1  e  Naroe : 

L2EG46 

r e 5  5 u r e  Con 

d  s  t i or  : 

V  a  c  u  u  m 

5  p  0 n  lv !  g  |  "•  C  1 

t  ; 

10.0  > m / s 

) 

a  t  a  Vw 

Tir 

Tout 

Ts 

£  <  rr,  •  £  > 

i  r  , 

(  C  ) 

(  C  ) 

1.70 

22.21 

26.12 

127.97 

2  1.7C 

22.21 

26.14 

128.12 

3  2.13 

22.03 

25.35 

127.99 

4  2.13 

22.08 

25.37 

128.06 

5  2.5! 

21  .99 

24.91 

128.07 

6  2.51 

21.99 

24.91 

123.00 

7  2. 89 

21.93 

24.57 

128.05 

6  2.89 

21.93 

24.56 

127.98 

9  3.32 

21.88 

24.25 

1 28 . 05 

0  4  .  J>2 

21.86 

24.25 

128.07 

t  3.55 

21  .85 

24.06 

127.99 

2  3 . 53 

21.85 

24.07 

128.03 

3  3.86 

21.82 

23.91 

128.03 

4  3.86 

21  .82 

23.90 

128.08 

in 

21.78 

23.64 

128.03 

6  4.40 

21.78 

23.65 

128.03 

7  1.70 

1  - - ? 

26.16 

128.03 

6  1.70 

<■*  ->  —  “ 
4.L  •  ki. 

26.15 

128.00 
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File 

Name  : 

L3EG5 1 

Pressure  Cor, a: 

t  ion : 

Vacuum 

Capo- 

Ce ioc 1 1  y 

: 

10.0  (1*1/5 

Data 

Uw 

Tin 

T  out 

Ts 

t 

( m/s  ) 

(C  1 

iC  ) 

(  C  ) 

1 

1  .70 

23.17 

27.31 

128.03 

-s 

1  .70 

23.18 

27.35 

128.03 

3 

2.13 

23.04 

28.52 

128.02 

4 

2.13 

23.05 

25.52 

128.05 

5 

2.51 

22 . 87 

25.05 

128.01 

6 

2.51 

22.97 

26.06 

128.08 

7 

2.88 

22.91 

25.73 

128.00 

S 

2.88 

22.91 

25.73 

127.99 

q 

3.31 

22.85 

25.40 

128.05 

10 

3.31 

22.88 

25.40 

128.03 

1  1 

3.58 

22.84 

25.24 

128.04 

1  2 

3 . 58 

22.84 

25.25 

128.00 

1  3 

3.85 

22.8! 

25.08 

128.04 

1  4 

3.85 

22.81 

25.08 

127.91 

15 

4.35 

'•v  *■>  •n  -> 

L.  .  {  ! 

24.82 

128.01 

15 

4. 35 

2 ~  7*7 

24.83 

128.01 

1  7 

1 .70 

23.2! 

27.38 

127.95 

1  £ 

1  .70 

23.2! 

n  ~~s 

127.57 

File 

Name : 

L4E850 

Pressure  C c 

ndit ion : 

Vacuum 

v  ape  r 

V  e  1  o  c 

‘  t  y  : 

i 0 . 0  (M'S 

1 

Data 

Cu 

7  in 

Tout 

Ts 

t 

'  m  s  / 

( C  :■ 

( C  ) 

( C  .• 

i 

1  .70 

23 . 00 

26 . 90 

128.02 

1  .70 

23.02 

26.97 

128.00 

J 

2.13 

22.90 

26.22 

128.05 

4 

2.13 

22.90 

26.20 

128.03 

c 

2.51 

22.82 

25.74 

127.97 

E 

2.51 

22.82 

25.75 

128.02 

7 

2.98 

22.77 

25.43 

128.01 

8 

2.88 

22.77 

25.43 

127.94 

9 

3.32 

22.72 

25.13 

128.00 

10 

3.32 

22.72 

25.13 

127.98 

1  1 

3.58 

22.70 

24.96 

127.93 

1  2 

3.58 

22.70 

24.97 

127.99 

13 

3.85 

22.63 

24.82 

128.02 

1  4 

3.85 

22.68 

24.83 

127.98 

15 

4.39 

22.64 

24.57 

127.97 

1  5 

4 . 39 

22.64 

24.57 

126.06 

17 

1  .70 

23.08 

27 , 08 

126.01 

1  S 

i  .70 

23.08 

"-.06 

128.02 
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Pile  Name:  L5E648 

Pressure  Condition:  Vacuum 
Vaoor  Velocity:  10.0  ( n/ s  ) 


Data 

Vw 

Tin 

Tout 

Ts 

t 

( m/  s  ) 

(C  ) 

(C) 

(C  ) 

1 

1  .70 

22.82 

26.55 

128.03 

n 

tm 

1  .70 

n "»  on 

<_  4.  .  OL 

26.55 

127.99 

3 

2.13 

22.69 

25.81 

127.93 

4 

2.13 

22.68 

25.82 

127.96 

5 

2.51 

22.60 

25.38 

127.98 

6 

2.5! 

22.60 

25.40 

128.01 

< 

2.89 

22.54 

25.04 

128.01 

8 

2 . 89 

22.54 

25.04 

128.00 

9 

3.32 

22.48 

24.74 

127.98 

10 

0.04 

22.49 

24.73 

127.96 

!  ! 

3.59 

22.46 

24.58 

128.01 

1Z 

3.59 

22.46 

24.58 

127.95 

13 

3.86 

22.43 

24.43 

127.97 

1  4 

3.86 

22.43 

24.44 

128.07 

15 

4.39 

22.35 

24.16 

128.00 

1  E 

4 .39 

22.35 

24.18 

127.94 

17 

1  .70 

22.64 

26.56 

127.94 

18 

1  .70 

22. S4 

26.54 

127.96 

Pile 

Marie : 

L6E&47 

Pres 

sure  Con 

di t l on : 

Vacuum 

Vapor  Ceioci 

♦  i .  • 

*■  J  * 

10.0  m  /  s 

> 

Data 

Vui 

Tin 

Tout 

T  5 

4 

\  M  /  S  • 

•,  c  > 

(C  :• 

<C) 

I 

1  .  70 

22.58 

25.90 

127.56 

2 

i  .  70 

22.56 

25.92 

125.03 

2.13 

22.46 

25.24 

127.92 

4 

2.13 

22.45 

25.24 

127.97 

c 

2.51 

22.39 

24.63 

127.99 

5 

2.51 

22.39 

24.81 

127.94 

7 

2.89 

22.34 

24.52 

127.95 

8 

2.89 

22.34 

24.54 

127.93 

9 

7  7  7 

■J  •  i. 

22.30 

24.26 

127.94 

10 

3.32 

22.30 

24.25 

127.98 

1  1 

3.59 

22.27 

24.05 

127. 96 

1  2 

3.55 

7  7  nn 

L.  fc.  .  L  i 

24.  10 

128.03 

1  3 

3.86 

o  C 
.  LO 

23.97 

127.96 

1  4 

3.86 

22 . 25 

23.97 

127.87 

15 

4 . 40 

0-1  0-3 

k.  •  4.  *. 

23.75 

127.96 

1  6 

4 . 40 

-in 

23.75 

127.93 

1  7 

1 .70 

22.67 

26.00 

127.94 

16 

1  .  70 

22.67 

25.01 

1Z7.SE 
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File  Name:  Q1RA20 

Pressure  Condition:  Atmospheric 
Uapor  Velocity:  .25  (m/s) 


Data 

Du 

Tin 

Tout 

Ts 

# 

(  m  /  s  ) 

(C  ) 

(C  ) 

(C  ) 

1 

2.07 

21  . 

61 

22.63 

48.52 

2 

2.07 

21  . 

61 

22.63 

48.54 

0 

2.64 

21  . 

43 

22.31 

48.60 

4 

2.64 

21  . 

43 

22.30 

48.56 

5 

3.51 

21  . 

26 

21.97 

48.51 

6 

3.51 

21  . 

26 

21.96 

46.54 

7 

4.47 

21  . 

15 

21  .73 

48.48 

8 

4 . 47 

21  . 

15 

21.73 

48.42 

9 

5.33 

21  . 

09 

21.58 

48.45 

)0 

5.33 

21  . 

09 

21.58 

48.41 

1  1 

6.  10 

21  . 

06 

21.49 

48.5! 

1  2 

6.10 

2  1  . 

06 

21.49 

48.46 

13 

6. 67 

21  . 

02 

21.40 

46.43 

1  4 

6.87 

2  1  . 

02 

21.4! 

48.4! 

15 

7.63 

20. 

96 

21.31 

48 . 35 

IE 

7 . 63 

20. 

96 

21.31 

46.39 

1  7 

2.07 

21  . 

75 

22.76 

46.74 

IS 

2 . 07 

21  . 

7C 

22.76 

46.77 

File 

Name : 

Q  STA106 

Pressure  Con 

di t ion : 

Atmospher ic 

Steam 

Ueioci 

ty  = 

1.0  ( m / s ) 

Date 

Dw 

T  l  n 

Tout 

Ts 

a 

(  m  '  s  ■ 

;  C  ; 

V  C  ) 

(C  ) 

1 

2.0E 

22.96 

30. BE 

100.03 

L. 

2.06 

22.96 

30.68 

99.98 

A 

•mJ 

2 . 64 

22.7/ 

29.48 

100.01 

4 

2.64 

"» n  nn 
l.  l.  .  <  r 

29.48 

100.05 

5 

3.31 

22.63 

28.55 

99.  S9 

6 

3.31 

22.63 

28.55 

100.02 

7 

3.98 

22.53 

27.85 

100.03 

8 

3.98 

22.54 

27.86 

100.0! 

9 

4.65 

22.46 

27.3? 

100.02 

10 

4.65 

22.46 

27.32 

100.06 

1  1 

5.32 

22.40 

26.88 

100.00 

1  2 

5.32 

22.40 

26.88 

100.00 

13 

5.99 

22.35 

26.50 

100.02 

14 

5.99 

22.35 

26.52 

100.02 

15 

6.47 

It.  . 

26.30 

100.00 

16 

6.47 

22 . 33 

26.30 

100.01 

17 

2.06 

23.05 

30.88 

100.03 

1  6 

2.06 

23.05 

30.87 

99.58 
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APPENDIX  B 


UNCERTAINTY  ANALYSIS 


When  measuring  a  physical  quantity,  there  will  always  be  a  differ¬ 
ence  between  its  actual  value  and  the  measured  value.  The  best  esti¬ 
mate  of  this  difference  is  defined  as  uncertainty  of  the  measured  value. 
The  uncertainty  depends  on  equipment  calibration  and  accuracy  as 
well  as  the  operator.  Although  an  uncertainty  for  a  single  measurement 
might  be  very  small,  an  equation  or  data  reduction  that  combines  two 
or  more  measurements  may  generate  results  with  rather  large  uncer¬ 
tainties.  The  uncertainties  of  the  physical  quantities  in  this  investiga¬ 
tion  are  combined  using  an  equation  suggested  by  Kline  and  Mclintok 
[Ref.  31]. 


Wr  = 


-,0.5 


(B.l) 


where 

Wr  =  uncertainty  of  the  desired  dependent  variable, 

xn  =  the  measured  variables,  and 

Wn  =  the  uncertainties  in  the  measured  variables. 

A  program,  designed  by  Mitrou  [Ref.  10]  to  compute  the  uncer¬ 
tainties  of  the  modified  Wilson  plot  technique,  was  used  in  this  inves¬ 
tigation.  Slight  modifications  were  made  to  include  the  difference 
inside  diameters  and  the  testing  of  ethylene  glycol. 


The  following  equation  was  used  to  compute  the  uncertainties  on 
the  the  enhancements. 


W 


c 


(B.2) 


where 

We  =  uncertainty  in  enhancement  for  finned  tube, 

Wrsm  =  uncertainty  of  outside  heat-transfer  coefficient  for  a 
smooth  tube  at  maximum  flow  rate, 

Wrsmn  =  uncertainty  of  outside  heat-transfer  coefficient  for  a 
smooth  tube  at  minimum  flow  rate, 

Wrfm  =  uncertainty  of  outside  heat-transfer  coefficient  for  a  finned 
tube  at  maximum  flow  rate,  and 

Wrfmn  =  uncertainty  of  outside  heat-transfer  coefficient  for  a  finned 
tube  at  minimum  flow  rate. 

The  following  are  examples  of  the  uncertainties. 
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DATA  FQR  THE  UNCERTAINTV  ANALYSIS: 


File  Name : 

Fnessure  Condition: 

S4RA17 

Atmospheri 

c  (101  i.  Pa  1 

Vapor  Temperature 

= 

48.42 

( Deg  C  ) 

water  Flow  Rate  <  X 

) 

= 

20.00 

Water  Velocity 

* 

2.08 

( m/s  > 

Heat  Flu-. 

6 . 333E+04 

(W/m'2 ) 

Tube-metal  thermal 

conduc . 

= 

385.0 

(U/m.k  ) 

Sieder-Tate  constant 

= 

0 . 0300 

UNCERTAINTY  ANALYSIS: 
VARIABLE 

Mass  F i ow  Rate,  lid 
Reynolds  Number,  Re 
Hea  t  Flu-,  q 
Log-Mean-Tem  D:ff,  LNTD 
Wail  Resistance,  R w 
0- era!  1  H . T . C .  ,  Uo 
Water-Side  H.T.C.,  Hi 
Uepor-Side  H.T.C.,  Ho 

DATA  FOR  THE  UNCERTAINTY 


PERCENT  UNCERTAINTY 

3.00 

3.12 
3.15 

.86 

2.63 

3.27 

7.13 
8.25 

ANALYSIS: 


File  Name :  S4RA 1 7 

Pressure  Condition:  Atmospher 

Ve.p o-  Tempe-ature  = 

Wate^  Flow  Rate  (  «  1  = 

wate-  Velocity  = 

Heat  Flu-  * 

Tube-metal  thermal  conduc.  * 

Sieder-Tate  constant  = 


C  v'10t  I  F  a  > 
48.40 
80 . 00 
7.83 

6 . 8S7E+04 
385.0 
0 . 0300 


(Deg  C  ) 

(  m  /  s  1 

(  V '  m  2  ! 
(  W  /  m  .  K  > 


UNCERTAINTY  ANALYSIS: 
VARIABLE 

Hass  F]0w  Rate,  Hd 
Reynolds  Number,  Re 
Heat  Flu.-  ,  q 
Log-Mean-Tem  Diff,  LHTD 
Wat  1  Resistance  ,  Rw 
Ove-a  1 1  H. T . C .  ,  Uo 
Wate^-Siae  H.T.C.,  Hi 
Vapor-Side  H.T.C.,  Ho 


PERCENT  UNCERTAINTY 

0.79 
1  .  IB 

3.13 
3.00 
2.63 
4.33 
6.74 
5.80 
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DATA  POR  Trt£  UNCERTAINTY  ANA_<£I£: 


File  Na*e: 

L4RA07 

Pressure  Ccndticr: 

AtnespHerjc  '.10!  tPa ) 

Vaoc-  Terce-atare 

-  46.43 

< Deo  C  ) 

Uate-  Flow  Rate  «  i  '■ 

-  20.80 

wate-  Velocity 

-  !  .16 

( n/s  ) 

Heat  Flu* 

•  8.07SE+04 

<U/W  2  ) 

TaBe-«etai  tAe-«al  conouc. 

-  385.0 

<  W/m.K  i 

Sieoer-Tate  constant 

-  0.0300 

UNCERTAIN* v  ANALYSIS: 

UARIABuE  PERCENT  UNCERTAINTY 

Mas?  F .Co  &ate  ,  A; 

2.00 

Re»rcies  Ny-te- .  ce 

3.10 

near  c *y •  .  z 

2.10 

L0s-Aear-Te* 

.65 

wall  ces  :  s  tarce  .  Rk- 

2.67 

C  .eral  1  jz 

3.16 

ware--: :  ce  **.  *  .C  .  .  ; 

'.accr-s •  se  h.  * .  2 .  .  h; 

i  7.46 

C“'~  =C=  *“£  UNCERTAINTY  ANAlySI 

5: 

F i .«  Ne-e:  »4RA0“ 

F-eit^-e  1  c  r  C .  r : :  -  :  -■  ; 

:  MC  t  Pa 

Vac:-  Te^se-et --e  » 

46.42  < Deo  C 

•are-  F;5fc  Fere  *  » 

60.00 

ware-  . e . : : ;  t .  • 

4.4C  *  * 

Hear  f.„.  * 

1.2' 2E-C5  w  -  2 

* e .  r -e- - * .  c:-.cy-.  » 

265 .  C  '  w  »  .f 

:;ece--*et e  :r-*ra-t  * 

0.02C0 

VN2 


vA»iABwE  PERCENT  UNCERTAINTY 


Hass  Flow  Fate ,  Ac 

0.79 

Feyrolas  Nunfcer .  Fe 

1.10 

Hea  t  f j u ,  t  a 

i  .  92 

l  03- Nee-- Ter.  Ds**,  LATD 

1.70 

Wall  Resistarce,  Rw 

2.6“ 

Overal I  m.T.c. ,  Vo 

2.5? 

Uate--Siae  h.T.C.  ,  hi 

6.72 

Vac  c--£i  ae  •*.  T  .C  .  , 

£.66 

137 


C'Ata  fQC,  ’hE  UNCERTA  I  NT  y  ANALYSIS  : 


File  Name:  Ml  ETC 103 

Pressu-e  Condition:  Vacuum  Ml  1  Pa  1 


Caper  Temperature  ■  46.44 

Water  Flow  Rate  '  X  >  *  20.00 

water  Velocity  ■  1.16 

Heat  Flu*.  -  1 . 352E+05 

Tude-metai  thermal  conduc.  ■  385.0 
Siede"-Tate  constant  *  0.0660 


<  DeQ  C  ) 

( m/ s  > 

(  W  /  m '  2  ) 
(  U! .  '  m .  K  ) 


UNCERTAINTY  ANALYSIS: 


CAP  I A2lE 


PERCENT  UNCERTAINTY 


Mass  FI  :w  Fate  ,  Md 

3.00 

Re>rcids  Number,  Pe 

2. 11 

neat  F* u>  ,  c 

3.06 

» c Mea r - T em  3 1  *  t  ,  lMT 

.41 

l  a  i  1  Resistance,  P  w 

2.E7 

C .  e-ai  i  r>.  * .  C  .  ,  L'c 

3.0? 

wete--£ : de  -. T . C  .  ,  Hi 

3.54 

v#rcr*: : ce  «. 7 . 3 .  ,  He 

on; 

Ca’A  F3c  Tr£  UNCERTA • NT » 

ANAL  t  S I E  : 

Fi le  Na-e:  M'S^VCc 

F'-essu-e  3c-dit.:-r:  Vacu-r  i 

i.ec:-  *e"ce"ftjre 
w»*e-  F  i c *  Rate  V 
Late'  -elicit. 

-eat  r  .  - ■ 

rube-*e  t  a  i  the-mei  eerduc  . 
E:eoe--*ate  carste-t 


UNlEP  * A I Nr •  ANAW*S1£: 

VARIABLE  PERCENT  UNCERTAINTY 


Mass  Flow  Rate  ,  Id 

0.96 

Reynolds  Number,  Re 

1  .24 

Hea  t  Flu*,  a 

1  .39 

Loa-Mear-Ter.  Lift,  LMTO 

.92 

Wail  Res i stance  ,  P* 

2.67 

Over e  1  i  w , T . 3 .  ,  Uo 

1  .67 

Water-Side  h.’.C.  ,  Hi 

3.21 

3.14 

>  >  Re  ■ 

48.44  >Ces  C 

E6 . 00 
3.64 

1 . B'CE’CE  1  W  n  3 
365 . 0 
0.0560 
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DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name '• 

S4STA100 

Fressure  Condition: 

Atmospher 

ic  (101  kPa  ) 

Vapor  Temperature 

= 

100.00 

( Deg  C  ) 

Water  Fiow  Rate  < % > 

20.00 

Water  Velocity 

« 

2.06 

( m/s  ) 

Heat  Flu* 

s 

5. 192E+05 

(W/i*T2  ) 

Tube-metal  thermal  conduc.  = 

385.0 

(  W  /  m .  K  ) 

Sieder-Tate  constant 

0.0510 

UNCERTAINTY  ANALYSIS: 
VARIABLE 


PERCENT  UNCERTAINTY 


Mass  Flow  Rate  ,  Md  3.01 
Reynolds  Number,  Re  3.14 
Heat  Flu-  ,  q  3.04 
Log-Mean-Tem  D l f  f  ,  LMTD  .11 
Wall  Resistance,  Rw  2.G3 
Overall  H.T.C.,  Uo  3.04 
Water-Side  H.T.C.  ,  Hi  6.41 
Vapor-Side  H.T.C. ,  Ho  2C.41 


DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name:  S45TA100 

Pressure  Condition:  Atmosphe'" 

Vapor  Temperature  * 

Wate-  Flow  Rate  • 

Water  velocity  = 

Hea t  F ] u .  = 

Tube-metal  thermal  corpuc.  *= 

Eiede--Tate  constant  = 


ic  <10!  Ira/ 
100.00 
66.00 
6.47 

7 . 693E+05 
38E.0 
0.0510 


(Deg  C  i 

i  m  /  s  ; 

(  W  '  m '  2  ) 
(  W  /  m  .  K  ) 


« 


UNCERTAINTY  ANALYSIS: 

VARIABLE 

Mass  Flow  Rate  ,  Md 
Reynolds  Number,  Re 
Heat  Flu*  ,  q 
uog-Meen-Tem  Diff,  LMTD 
Wal 1  Resistance  ,  Rw 
Over  a i  1  H.T.C.,  Uo 
Water-Side  H.T.C.,  Hi 
Uapor-Side  H.T.C.,  Ho 


PERCENT  UNCERTAINTY 

0.96 
1  .31 

1  .06 
n  n 

2.63 

1.10 

5.99 

6.51 
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DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name: 

SMTHSTA65 

Pressure  Condition: 

Atmospheric  (101  KPa ) 

Vapor  Temperature 

= 

99.96 

( Deg  C  ) 

Water  Flow  Rate  <  %  ) 

= 

20.00 

Water  Velocity 

s 

1.16 

( m/  s  ) 

Heat  Flux 

= 

4 . 003E+05 

(  U  /  m  '■  2  ) 

Tube-metal  thermal  conduc. 

= 

385.0 

( W/m.K  ) 

Sieder-Tate  constant 

= 

0.0660 

UNCERTAINTY  ANALYSIS: 


VARIABLE 


PERCENT  UNCERTAINTY 


Mass  Flow  Rate,  Md  3.00 
Reynolds  Number,  Re  3.11 
Heat  Flu*,  q  3  .'04 
Log-Mean-Tem  Diff,  LMTD  .14 
Wall  Resistance,  Rw  2.G7 
Overall  H.T.C.,  Uo  3.04 
Water-Side  H.T.C.,  Hi  3.94 
Vapor-Side  H.T.C.,  Ho  7.91 


DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name: 

Pressure  Condition: 
Vapor  Temperature 
Water  Fiow  Rate  ' %  > 
Water  Velocity 
meat  Flu-. 


SMTHSTAB5 

Atmospheric  <10!  |. Pa  > 
=  100.02 
=  80.00 
=  4.40 

=  5.301E+05 


Tude-netal  thermal  conduc.  =  385.0 

Sieder-Tate  constant  =  0.0E50 


<  Deg  C  ! 

( m  /  s  i 
(  W  /  m  ’  2  ) 

<  W  /  m .  K  ) 


UNCERTAINTY  ANALYSIS: 
VARIABLE 

Mass  Flow  Rate,  Md 
Reynolds  Number,  Re 
Heat  Flu*  ,  q 
Log-Mean-Tem  Diff,  LMTD 
Wall  Resistance  ,  Rw 
Overa  1 1  H . T . C .  ,  Uo 
Water-Side  H.T.C.,  Hi 
Vapor-Side  H.T.C.,  Ho 


PERCENT  UNCERTAINTY 

0.79 
1.12 
.99 
.39 
2.67 
1  .06 
3.18 
1.93 
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DA,TA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name:  L4STAE6 

Pressure  Condition:  Atmospher 

1C  (101  kPa ) 

Vapor  Temperature 

c 

99.95 

(Deg  C) 

Water  Flow  Rate  (%) 

= 

20.00 

Water  Velocity 

- 

1.16 

(  m  /  s  > 

Heat  Flux 

= 

6.397E+05 

<  W/mA2 ) 

Tube-metal  thermal  conduc. 

s 

385.0. 

( W/m . K  ) 

Sieder — Tate  constant 

SB 

0.06E0 

UNCERTAINTY  ANALYSIS: 


VARIABLE 


PERCENT  UNCERTAINTY 


Mass  FIooj  Rate,  Md  3.00 
Reynolds  Number,  Re  3.12 
Heat  Flux,  q  3. -04 
Log-Mean-Tem  Diff,  LMTD  .09 
Wall  Resistance,  Rw  2.B7 
Overall  H.T.C.,  Uo  3.04 
L'ater-Side  H.T.C.,  Hi  3.94 
Vapcr-Side  H.T.C.,  Ho  25.55 


DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name:  L4STAEE 

Pressure  Condition:  Atmospheric  (101  kPa  ) 

Vapor  Temperature  =  100.01 

Later  Flow  Rete  (T>  *  80.00 

Water  Velocity  =  4.40 

Heat  Flux  =  1.122E+06 

Tube-metal  thermal  conduc.  =  385.0 

Sieder-Tate  constant  =  0.0560 


(Deg  C> 

<1*1/5  ) 

( W  /  m "  2  ) 
(  W  /  m .  K  5 


< 


\ 


UNCERTAINTY  ANALYSIS: 

VARIABLt  PERCENT  UNCERTAINTY 


Mass  Flow  Rate  ,  Md 

0.79 

Reynolds  Number,  Re 

1.12 

Heat  Flux,  q 

.93 

Log-Mean-Tem  Diff,  LMTD 

.19 

Wall  Resistance,  Rw 

2.67 

Overall  H.T.C.  ,  Uo 

.94 

Water — Side  H . T . C .  ,  Hi 

3.18 

Vapor — Side  H .  T . C .  ,  Ho 

5.18 
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DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name: 

L4STVE7 

Pressure  Condition: 

Vacuum  < 

1  1 

t  Pa  ) 

Vapor  Temperature 

= 

48.49 

( Deg  C  > 

Water  Flow  Rate  ( %  ) 

- 

20.00 

Water  Velocity 

= 

1  .  i  6 

(  m  /  s  ) 

Heat  Flux 

•=c 

1 . 885E+05 

(  W  /  m  A  2  ) 

Tube-metal  thermal  conduc. 

385.0 

(W/m.K  ) 

Sieder-Tate  constant 

= 

0.0EG0 

UNCERTAINTY  ANALYSIS: 


VARIABLE 

PERCENT  UNCERTAINTY 

Mass  Flow  Rate,  Md 

3.0! 

Reynolds  Number  ,  Re 

3.  I  ! 

Heat  F lu.*.  ,  q 

3. '05 

Log-Mean-Tem  Diff,  LMTD 

.30 

Wall  Resistance,  Rw 

2 .87 

Over a  1 1  H . T . C .  ,  Uo 

3.07 

Water — Side  H.T.C.  ,  Hi 

3.94 

Vapor-Side  H.T.C.  ,  Ho 

23. 9S 

DATA  FOR  THE  UNCERTAINTY 

ANALYSIS: 

File  Name:  L4STVB7 

Pressure  Condition:  Vac 

uum  (  1  1 

kPa  > 

Vapor  Temperature 

= 

48.49 

<  Deg  C 

W  a  t  e  -  Flow  Rate  <  %  ■ 

s 

80 . 00 

Water  Velocity 

a 

4.35 

(  mi  /  s  ' 

Heat  Flux 

* 

3 . 779E+05 

(  W  /  m '  Z 

Tube-metal  thermal  conduc 

* 

385.0 

(  W  /  m  .  F 

Sieder-Tate  cons t ant 

= 

0 . 0580 

UNCERTAINTY  ANALYSIS: 

VARIABLE  PERCENT  UNCERTAINTY 


Mass  Flow  Rate  ,  Md 

0.79 

Reynolds  Number,  Re 

1.12 

Heat  Flux,  q 

1  .08 

Log-Mean-Tem  Diff,  LMTD 

.55 

Wall  Resistance  ,  Rw 

2.87 

Overal 1  H.T.C.,  Uo 

1.19 

Water-Side  H.T.C.  ,  Hi 

3.18 

Vapor — Side  H.T.C.,  Ho 

8.58 
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HAT  A  CHE  Ti-!r  f  iMHrE  T  A  T  |\{T  v  AWA!  VCTQ* 


MPPGU74 


Pressure  Ccn  d :  t  i  or*  r  u=c  (  *  -  VPs) 


Uar.r.r  Taiv»nar‘a  +  . 


__  in'?  tq 


( n»r.  r  > 


L!a  +  ar  C"1  DaU  ^ 


=  50.05 


L*s+a:r  W«  1  .n-  • ■  +  . 


+  PI, 


C  i 1 1 C+flC  /  !J  /  iv.  *-  n  v. 


T  uL1  5  ~  W  £  t  5  1  t  K  a r~  ivs  a  1  C  C*  f*  d  l 


=  335.0 
=  0.0350 


ij  /  m  1 f.  ) 


!  INf'FRTATMTY  AMAIS'FTF: 


(JAR  JABLE 


PERCENT  UNCERTAINTY 


Ns;;  Flow  Rate,  tfd 

P.& v,.  r%.~.  i  h ~  Number  Re 


_ _MC^,._T.W  p|  ,  f  f  l  mt n 


!.!=•’  5ic:=»;r,.-;  Di., 


n..Bra1  !  U  T 


Nste^-Side  U.T.r^  Hi 


Uar.-,.--R  •  rt* 


qata  FQP  the  UNCERTAINTY  ANALYSIS? 


F-.  !*  Name 


MFF£U40 


qata  fqp  tuc  UNCERTAINTY  ANALYSIS: 


F  •  i  a 


MFFOU~4 


Prejt.jre  Condi  t  i  0 r- •'  U5C  UUfB  •'  1  >  !  pi  i 


=  RC  {*!» 


I,!  5  *  I.'a  ’  •,  t  . 


PI,,. 


f  RFAF *r»r 


T  ;  ,  t  -  ,v  ^  i  z,  1  +  t .  £5  r-  r,*,  a  '  r  r.r,  H .  .  r 


t-.eHer-Tsle  fr.r.ciaf 


UNCERTAINTv  ANALYSIS? 


war  table 


PERCENT  UNCERTAINTY 


Nese  Flow  Rate,  M  d 
R e y no Ida  N u m fc e r  ,  Re 
Heat  Flu,,  q 
LoQ-Plean-Tem  Diff,  LP1TD 
Wall  Resistance ,  R w 
Qverg] 1  h.T.C. ,  Uo 
w*+  — P i  rt»  m  t  r  u> 

-  -  —  ■  ~  ~  i 

vaqc‘--Side  H.T,C.  Me 


DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name: 

M1EGV37 

Pressure  Condition: 

Vacuum 

(  1  1 

IP  a  ) 

Vapor  Temperature 

= 

127.83 

( Deg  C  ) 

Water  Flow  Rate  ('/.) 

= 

20.00 

Water  Velocity 

= 

1  .  18 

<  m/ s  ) 

Heat  Flux 

= 

2. 1 6BE+05 

<W/m"2  ) 

Tube-metal  thermal  conduc. 

= 

385.0 

<  W/m . K  ) 

Sieder — Tate  constant 

= 

0.0350 

UNCERTAINTY  ANALYSIS: 


VARIABLE 


PERCENT  UNCERTAINTY 


Mass  Flow  Rate,  Md  3.00 
Reynolds  Number,  Re  3.11 
Heat  Flux  ,  q  3  .'04 
Log-Mean-Tem  Diff,  LMTD  .25 
Wall  Resistance,  Rw  2.6? 
Ov'eraii  H.T.C.,  Uc  3.05 
Uater-Siae  H.T.C.,  Hi  14.51 
Vapor-Si de  H. T . C . ,  Ho  10.10 


DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name:  M1EGV37 

Pressure  Condition:  Vacuum  > 

Vapor  Tempe-ature 
Later  Flow  Rate  < %  ) 

Water  Velocity 
Heat  Flux 

Tube-metal  thermal  conduc. 
Sieder-Tate  constant 


KPa  > 

127.81  'Deg  C : 

80 . 00 

4.40  ( m/s  > 

b . 50 1 E^CB  ( W / m ’ 2  ! 
385 .0  (  W m  .  K  / 

0 . 0350 


UNCERTAINTY  ANALYSIS: 
VARIABLE 

Mass  Flow  Rate,  Md 
Reynolds  Number,  Re 
Heat  Flux  ,  q 
Log-Mean-Tem  Diff,  LMTD 
Wall  Resistance,  Rw 
Overal 1  H . T . C .  ,  Uo 
Water-Side  H.T.C.,  Hi 
Vapor-Side  H.T.C.,  Ho 


PERCENT  UNCERTAINTY 

0.79 
1.11 
1  .22 
.82 
2.67 
1  .48 
14.32 
3.39 
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OAT A  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name:  S4STU74 

Pressure  Condition:  Vacuum  (1 

Uapor  Temperature  = 

Wa t er  F 1 ow  Rate  ( %  )  = 

Water  Velocity  = 

Heat  Flux  = 

lube-metal  thermal  conduc.  = 

Sieder-Tate  constant  = 


KPa  ) 

48.45 

20.00 

2.06 

1 .592E+05 
385.0 
0.0510 


(Deg  C) 

<  m/s  ) 

( W/rTZ  > 
(W/m.K  ) 


UNCERTAINTY  ANALYSIS: 


UARI ABLE 

Mass  Flow  Rate,  Md 
Reynolds  Number,  Re 
Heat  Flux  ,  q 
Log-Mean-Tem  Oiff,  LMTD 
Wall  Resistance,  Rw 
Overall  H.T.C. ,  Uo 
Water-Side  H.T.C.  ,  Hi 
Uapor-Side  H.T.C.,  Ho 


PERCENT  UNCERTAINTY 

3.00 
3.  13 
3.05 
.35 
2.63 
3.07 
6.41 
18.78 


DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name:  S4STU74 

Pressure  Condition:  Uacuum  (1 

Vapor  Temperature  = 

Water  Flow  Rate  (%)  = 

Water  Velocity  = 

Heat  Flux  * 

Tube-metal  thermal  conduc.  = 

Sieder-Tate  constant  = 


KPa) 

48. 4G 
40 .00 
3.98 

2.01 IE+05 
385.0 
0.0510 


(Deg  C) 

( m/s  ) 

( W/n'2 ) 
(W/m.K  ) 


UNCERTAINTY  ANALYSIS: 
VARIABLE 

Mass  Flow  Rate,  Md 
Reynolds  Number,  Re 
Heat  Flux,  q 
Log-Mean-Tem  Diff,  LMTO 
Wall  Resistance,  Rw 
Overall  H.T.C.  ,  Uo 
Water-Side  H.T.C.,  Hi 
Vapor-Side  H.T.C.,  Ho 


PERCENT  UNCERTAINTY 

1  .55 
1  .78 
1  .70 
.53 
2 .63 
1 .78 
6.07 
8.85 
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DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name :  RSMTH1 

Pressure  Condition:  Atmospheric 

Vapor  Temperature  * 

Water  Flow  Rate  (X)  = 

Water  Velocity  = 

Heat  Flux  = 

lube-metal  thermal  conduc.  =  3 
Sieder-Tate  constant  = 


( 101  IPa ) 
49.22 
20.00 
1 .17 

2.S91E+04 
85.0 
0 . 0300 


( Deg  C  ) 

<  m/s  ) 

( W/m',2  ) 
(U/m.K  ) 


UNCERTAINTY  ANALYSIS: 
VARIABLE 

Mass  Flow  Rate  ,  Md 
Reynolds  Number,  Re 
Heat  Flux,  q 
Log-Mean-Tem  Diff,  LMTD 
Wall  Resistance,  Rui 
Overall  H.T.C.,  Uo 
Water-Side  H.T.C.,  Hi 
Vapor-Side  H.T.C.,  Ho 


PERCENT  UNCERTAINTY 

2.39 
3.08 
3. 58 
2.11 
2. 57 
4.24 
7.12 
5.71 


DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name:  RSMTH1 

Pressure  Condition:  Atmospheri 

Vapor  Temperature  * 

Water  Flow  Rate  <%) 

Water  Velocity  = 

Heat  Flux  = 

Tube-metal  thermal  conduc.  = 

Sieder-Tate  constant  = 


c  <101  1  Pa  ) 
48.45 
80 . 00 
4.42 

2 . 995E+04 
305.0 
0 . 0300 


(Deg  C) 

<  m  /  s  ) 
(W/m-'Z) 
(U/m.K  > 


UNCERTAINTY  ANALYSIS: 
VARIABLE 

Mass  Flow  Rate,  Md 
Reynolds  Number,  Re 
Heat  Flux  ,  q 
Log-Mean-Tem  Diff,  LMTD 
Wall  Resistance,  Ru 
Overall  H.T.C.  ,  Uo 
Watei — Side  H.T.C.  ,  Hi 
Vapor-Side  H.T.C.,  Ho 


PERCENT  UNCERTAINTY 

0.79 

1.08 

6.98 

6.92 

2.67 

9.82 

6.73 

10.83 
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